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Abstract 

The  statistical  mechanics  of  stationary,  homogeneous,  hydromaRnetic 
turbulence  is  treated  by  a  method  which  does  not  require  the  assumption  of  small 
departure  from  a  state  of  detailed  balance.  The  equations  of  motion  under  cyclic 
boundary  conditions  are  expressed  in  terms  of  the  external  driving  forces  and 
linearly  independent  dynardc  variables  formed  from  the  spatial  Fourier  coefficients 
of  velocity  and  magnetic  fields.  The  nonlinear  interaction  is  exhibited  as  the  sum 
of  individually  conservative  elementary  interactions  among  triads  of  modes.  A 
functional  equation  is  constructed  giving  necessary  and  sufficient  conditions  for 
a  distribution  of  time  functions  to  satisfy  the  equations  of  motionj  it  involves 
only  the  second-,  third-,  and  fouorth-order  moments  of  the  distribution  functional. 
By  using  the  fact  that  each  mode  interacts  with  manj'  others,  a  procedure  is 
developed  in  which  each  elanentsry  interaction  is  introduced  as  a  perturbation  on 
the  total  motion.  The  moments  appearing  in  the  functional  equation  are  expressed 
by  means  of  this  procedure  in  terms  of  the  diagonal  elements  of  the  covariance 
matrices  of  the  dynamic  variables  and  external  forces  and  the  distribution- 
averaged  infinitesimal  impulse  response'  matrix  of  the  system.  Closed  equations 
for  these  quantities  are  obtained  so  that  the  power  spectra  in  frequency  and 
wave  number  of  velocity  and  magnetic  fields,  and  the  energy  exchange  anong  the 
modes,  may  be  determined  from  the  driving  force  power-spectrum.  A  variational 
method  is  given  for  defining  and  constructing  the  most  nearly  normal  distribution 
functional  with  the  previously  derived  second-,  third-,  and  fourth-order  moments, 
thereby  determining  the  higher  moments.  The  consistency  with  the  present  theory 
is  confirmed  for  the  cascade  hypothesis  fundamental  to  the  Kolmogorov  sijnilarity 
theory  of  the  inertial  range.  It  is  found  that  the  rate  of  direct  energy-transfer 
from  modes  of  wave  number  below  k,  to  modes  of  wave  number  above  k^,  »  k,  is 
asymptotically  proportional  to  e(k, /k^)  '•^,  where  e  is  the  total  rate  of  energy 
dissipation  per  unit  mass  and  k,  ,kp  are  both  within  the  inertial  range. 
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Corrections 

The  normalization  of  the  integral  transform  U(kjr)  introduced  in  (9.8)  is 
correctly  given  by 

U(k,r)  =  (2n)"3  J  U(x,r)  e'^^*2^d^  X  . 
In  (A»ll)>  instead  of  -j()  read  +j()  • 

Additional  Comments 

The  corroboration  of  the  Kolmogorov  similarity  theory  given  in  Section  9  has 
a  doubtful  significance!  It  should  be  noted  that  the  universal  form  for  the  r  • 
dependence  of  U(k, r)  given  in  (9*12)  is  physically  a  rather  implausible  extension 
of  dimensional  reasoning.  The  convective  action  of  the  large  eddies  may  be  ex- 
pected to  affect  the  time-dependence  at  high  wave  numbers ^  so  that  even  if  U(k,0) 
actually  is  a  universal  function  in  the  equilibrium  range,  unaffected  by  direct 
interaction  with  the  large  eddies,  U(k,  t)  probably  is  not.  (This  effect  of  con- 
vection has  been  discussed  by  Kolmogorov  |^2]  and  ■  Heisenberg  [lO]«)  Although  (9*12) 
yields  a  consistent  viniversal  expression  for  S(k,0),  as  indicated  in  the  text,  it 
does  not  follow  from  this  that  the  expression  for  S(k,  r),  ("^Z  O),  is  independent 
of  the  low-lying  irave  numbers.  Thus,  the  consistency  of  (9»12)  vri.th  the  full 
equations  of  the  present  theory  has  not  been  demonstrated.  "" 

Further  consideration  of  the  convective  effects  of  the  ]arge  eddies  leads  to 
the  conclusion  that  the  idealized  model  of  the  present  paper,  in  which  the  enel-gv- 
containing  range  is  maintained  by  the  action  of  homogeneous  driving  forces  F(k,t) 
that  are  statistically  independent  for  different  k's,  may  differ  from  laboratory 
turbulence  more  significantly  than  was  indicated  in  Section  10.  In  the  model  the 
dimensions  of  the  system  are  very  large  compared  to  the  reciprocal  of  a  typical 
wave  number  in  the  principal  energy-containing  region,  and  the  phases  at  adjacent 
wave  numbers  in  this  region  are  essentially  unrelated  in  a  typical  realization  of 
the  process.  In  laboratory  turbulence  the  localization  of  the  system  to  dimensions 
comparable  to  eddy  sizes  in  the  energy-containing  region  (smd/or  the  generation  of 
the  large  eddies  by  fixed  mechanical  structures)  indicates  that  the  phase  in  this 
region  varies  only  smoothly  with  wave  number  in  a  typical  realization.  Thus  the 
low-lying  Fourier  modes  in  laboratory  turbulence  exhibit  a  snecies  of  coherence 
^■ihich  is  lacking  in  the  model.  The  presence  or  absence  of  this  coherence  conceiy^ 
ably  may  affect  interactions  at  the  higher  wave  numbers- by  influencing  the  extent 
to  wlxLch  the  low-lying  Fou^'-.r-r  modes  effectively  convect  the  high  wave  number 
region  as  a  whole,  or  disturb  internal  phase  relations  in  that  region. 

June  1957 
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1,  Introduction 

At  the  Reynolds  numbers  characteristic  of  turbulence  there  exist,  for 
given  boundary  conditions,  a  rich  manifold  of  distinct  solutions  of  the  hydro- 
dynamic  equations  of  incompressible  motion'-  -^ .  The  individual  solutions  have 
very  complicated  structure  and  are  unstable,  going  over  into  each  other  under 
the  action  of  perturbations.  We  often  are  interested  in  properties  of  the  motion 
which  can  be  expressed  as  averages  over  times  long  compared  with  characteristic 
periods  of  the  individual  solutions.  In  view  of  the  presence  in  any  physical 
situation  of  irreducible  random  perturbations  of  the  flow  it  is  reasonable  to 
assume  that  such  properties  can  be  identified  with  statistical  mechanical 
averages  over  an  appropriate  distribution  of  individual  solutions. 

The  real  and  imaginary  parts  of  the  spatial  Fouirler  coefficients  of  the 
velocity  field  (and  magnetic  field,  if  present)  may  be  regarded  as  phase-space 
coordinates  analogous  to  the  canonical  coordinates  and  momenta  of  Hamiltonian 
systems.  It  then  may  be  shown  that  the  system  obeys  a  Liouville  equation  similar 
to  that  for  Hamiltonian  systems'-  -' .  Ther«  are  essential  differences,  however, 
between  a  turbulent  fluid  and  systems,  such  as  a  dilute  gas  of  hard  spheres,  for 
which  the  Gibbs  statistical  mechanics  has  been  most  successful. 

In  a  dilute  gas  of  hard  spheres,  the  particles  are  noninteracting  almost 
all  of  the  time,  and  the  interaction  can  be  well  described  in  terms  of  transi- 
tion probabilities,  Hue  to  two-particle  collisions,  between  states  of  noninter- 
action. The  collisions  represent  a  small  time-averaged  perturbation  in  the 
Hamiltonian  of  the  uncoupled  particles.  In  the  turbulent  fluid,  a  tjrpical 
Fourier  mode  undergoes  simultaneous  and  continuous  interaction  with  a  large 
number  of  other  modes.  The  motion  is  due  entirely  to  the  interaction.  Fbr  the 
dilute  gas  we  may  be  interested  principally  in  the  velocity  distribution  of  the 
particles  and  interested  in  the  collisions  only  in  terms  of  their  frequency  and 
scattering  angles.  For  the  fluid,  since  the  interaction  is  continuous,  we 
necessarily  are  interested  in  the  behavior  of  the  Fourier  modes  during  the  inter- 
action. The  detailed  structure  of  the  interaction  cannot  be  ignored. 

We  may  note  also  that  stationary  turbulence  represents  a  very  strong  depar- 
ture from  a  state  of  detailed  balance.  It  follows  from  Kolmogorov's  similarity 

[2!* 
considerations'-  J  that  the  temperature  characterising  the  distribution  of  energy 

among  the  individual  modes  in  the  inertial  range  varies  with  wave  number  as  k~  '-^. 

Thus  it  cannot  be  expected  that  methods  of  irreversible  statistical  mechanics 

*See  also  G.K,  Batchelor^ -J . 
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based  on  small  departure  from  detailed  balance  could  be  successful  in  describing 
the  transport  of  energy  across  the  wave  n-omber  spectrum, 

A  third  difference  between  a  turbulent  fluid  and  a  dilute  gas  is  that  the 
Fourier  modes  of  the  fluid  all  enter  differently  into  the  equations  of  motion. 
The  distribution  function  must  exhibit  a  much  smaller  degree  of  degeneracy  than 
that  for  a  gas  of  identical  particles. 

G.  I.  Taylor'-  -'j  followed  by  many  other  authors'--^,  has  formulated  the 
statistical  problem  of  turbulence  as  that  of  determining  moments  of  the  velocity 
distribution.  The  moments  of  greatest  interest  are  contained  in  the  space-time 
covariance  tensor  of  the  velocity?  this  determines  the  spectral  distribution  of 
the  kinetic  energy  in  wave  number  and  frequency.  Equations  of  motion  for  various 
distribution  moments  result  from  multiplying  the  Navier^tokes  equation  with 
suitable  functions  and  averaging.  As  a  consequence  of  nonlineaidty,  the  equation 
for  the  velocity  covariance  contains  the  third  moments  also,  and  direct  attempts  at 
elimination  of  the  latter  lead  to  a  series  of  incomplete  systems  of  equations 
Involving  successively  higher  moments, 

A  part  of  this  grave  difficulty  seems  due  to  incomplete  formrilation. 
For  it  is  not  to  be  expected  that  the  Navier-Stokes  equation  determines  a  unique 
form  of  statistical  distribution,  even  under  the  conditions  of  isotropy  and 
stationarity,  unless  one  specifies,  directly  or  indirectly,  the  full  statistical 
structure  of  the  forces  irtiich  supply  energy  to  the  turbulence.  To  make  this 
plausible,  we  may  note  that  in  a  conservative  Hamiltonian  system  the  specification 
of  mean  energy,  and  the  imposition  of  stationarity,  leave  the  higher  momenti 
largely  undetermined,  since  the  simultaneous  distribution  may  be  any  function  of 
the  constants  of  motion.  If  dissipation  and  driving  forces  are  introduced  the 
ambiguity  still  remains,  unless  the  full  statistical  distribution  of  the  driving 
forces  is  fixed. 

Whatever  supplementary  conditions  are  imposed  on  the  distribution  must 
be  consistent  with  the  equations  of  motion  and  appropriate  to  the  physical  situa- 
tion to  be  represented, 

2,  Description  of  the  present  approach 

The  present  paper  developed  from  the  following  argument,  which,  for 
brevity,  we  state  for  nonmagnetic  turbulence.  Let  the  velocity  field  u,  be 
represented  in  the  usual  manner  by  Fourier  modes.  The  interaction  among  three 
modes  u(k),  u(k'),  u(k-k ' )  can  proceed  slon^   a  direct  path,  represented  by  terms 
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involving  u(k' )u(k-k' )>  u(k)u(k'-k),  u(k)u(-k')  in  the  respective  equations  of 
motion  of  these  modes,  and  along  indirect  paths  involving  chains  of  interactions 
with  the  other  excited  modes.  In  view  of  the  mixing  action  characteristic  of 
turbulence  let  us  suppose  that  the  combined  action  of  all  the  indirect  inter- 
actions does  not  convey  phase  information  among  the  three  modes,  so  that  the 
triple  correlation  among  them  is  due  to  the  presence  of  the  direct  interaction. 
A  pertinent  analogy  would  be  the  mutual  phase  relations  induced  by  the  interaction 
of  three  nonlinear  oscillators.  If  many  modes  are  excited  and  interacting,  the 
direct  interaction  shovild  represent  a  small  disturbance  of  the  motion  of  the  three 
modes,  and  then  it  should  be  possible  to  calculate  by  perturbation  methods  the 
triple  correlations  which  appear  in  the  equation  of  motion  of  the  spectrum  tensor. 

In  the  following  sections  it  is  shown  that  a  self-consistent  procedure  of 
this  sort  can  be  carried  out  in  which  each  elementary  interaction  among  three 
Fourier  modes  is  introduced  as  a  perturbation  on  the  total  motion  due  to  all  the 
interactions.  It  results  in  closed  equations  which  fix  the  covariance  matrix  of 
the  Fourier  modes  and  a  matrix,  basic  to  the  perturbation  treatment,  describing 
the  impulse  response  of  the  distribution  to  small  disturbances. 

The  statistical  mechanical  method  used  involves  a  distribution  functional 
of  time  functions  satisfying  the  equations  of  motion.  It  is,  in  this  respect, 
more  general  than  the  Gibbs  statistical  mechanics  ivhich  treats  a  distribution  only 
of  simultaneous  amplitudes,  and  which,  as  we  shall  see,  is  inadequate  for  the 
stationary  turbulence  problem.  Apart  from  its  specific  usefulness  in  the  turbulence 
problem,  the  method  presented  is  of  interest  because  its  applicability  is  not 
limited  by  the  requirement  that  the  state  of  the  system  be  close  to  a  state  of 
detailed  balance. 

The  ambiguity  in  the  distribution  functional  mentioned  in  the  Introduction 
is  resolved  by  the  requirement  that  the  distribution  of  the  Fburier  modes  be  as 
independent  as  the  equations  of  motion  permit.  This  is  expressed  first  implicitly 
by  a  hjrpothesis  of  weak  statistical  dependence  of  the  modes  (in  a  sense  specifically 
defined)  and  later  explicitly  by  a  variational  method  of  determining  the  higher 
moments  of  the  distirLbution  functional, 

3.  Representation  of  the  fields  by  linearly  independent  variables 

The  equations  of  motion  of  an  incompressible,  highly  conducting,  and 
noni^lativistic  hydromagnetic  fluid  subjected  to  external  body  or  boundary  forces 
may  be  written '--'J 


-  u- 

(3.1)  \{k)  +  vk^u^(k)  -  ik^  P^j(k)  2!  k»   [*3(iS-]S')wj^  (Is')  -  ^j(iS-iS')uj^  (iS')l 

*  Pi/k)  Fj(k), 

(3.2)  w^(k)  +  vk^w^(k)  -  ikj  2!kt    K^iS-lS'^Wj^iS')  -  '^j(iS-iS'K(k')l. 

Here  u(k),  w(k),  F(k)  ar«  related  to  the  velocity  field  u(x),  magnetic  iiiduction 
B(x),   and  external  force  field  F(x)  by 

(3.3)  u(x)  =  T  u(k)e^,     B(x)  =  {hnup)^^^  T  Me^,  F(x)  -pT  F(k)e^, 

and  P.  .(k)  ;  6 .  .  -  k   k.  k..  The  summations  range  over  all  wave  vectors  allowed 

ij  ""  ~  j-j       ^  J  _        -12 

by  the  boundary  conditions.  The  quantities  p,  v,  |i,(r,  c  and  v  »  (Uk|j,o~)   c 

denote,  respectively,  density,  kinematic  viscosity,  permeability,  conductivity, 

velocity  of  light  in  vacuo,  and  'ohmic  viscosity'}  they  all  are  assumed  constant 

throughout  the  fluid.  The  fields  also  obey  the  divergence  conditions 

(3.U)      k^  u^  (k)  -  0,    k^  w^  (k)  -  0, 

which  are  preserved  by  (3.1) »  (3.2),  if  imposed  initially. 
The  total  energy  per  unit  mass  is 


(3.5)     ^  "  I  Zk  K^IS)  ^i(k)  *  wj(^)  w^(k)j. 


Energy  is  conserved  by  the  nonlinear  interaction  of  the  Fourier  modes  in  the 
sense  that  if  v,  v  and  F(^)  are  put  equal  to  zero,  (3.l)»  (3.2),  (3.'4)  yield 
E  =  0. 

We  shall  assume  that  the  velocity  and  magnetic  fields  obey  cyclic  boundary 
conditions  on  the  faces  of  a  cube  of  side  L,  These  conditions  are  more  natural 
for  the  description  of  homogeneous  statistical  distributions  than  the  requirement 
that  the  fields  vanish  on  the  boundaries,  and  the  small-scale  dynamics  should  be 
insensitive  to  the  precise  form  of  the  boundary  conditions.  The  allowed  k-vectors 
now  are  all  those  whose  components  along  the  coordinate  axes  are  integer  multiples 
of  2n/L.  We  shall  further  assvime,  however,  that  there  is  no  uniform  velocity  or 
magnetic  field,  so  that  u(0)  =  w(0)  =0.  These  conditions  are  preserved  by  the 
equations  of  motion,  provided  j;(0 )  =  C. 
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The  coefficients  u(k),  w(k)  for  allowed  k  are  not  linearly  independent; 
assignment  of  their  initial  values  is  restricted  by  (3.U)  and  the  reality  require- 
ments 

(3.6)  u^(-k)  =  u*(k),      w^(-k)  -  wj(k). 

In  the  statistical  mechanical  treatment  it  is  a  great  advantage  to  woi^c  with 
variables  for  which  an  arbitrary  assignment  of  initial  values  corresponds  to 
admissible  fields.  For  this  purpose  we  introduce  those  real  and  imaginary  parts 
of  the  vector  components  of  the  u(k),  w(k)  which  are  linearly  independent  under 
(3.U),  (3.6). 

The  vector  components  of  u(k),  u(-^),  w(k),  w(-]5)  comprise,  for  given Jc, 
a  total  of  twenty-four  real  and  imaginary  parts,  of  which  only  eight  are  linearly 
independent.  We  shall  pick  from  each  wave  vector  pairjk,  ~k   one  member,  say  k, 
and  choose  as  these  eight  variables  the  real  and  imaginary  parts  of  the  components 
of  the  vectors  u(k),  w(k)  along  pairs  of  perpendicular  axes  in  the  plane  normal 
to  k.  Let  all  the  independent  real  and  imaginary  parts  so  chosen  (for  all  pairs 
k,  -^   allowed  by  the  boundary  conditions)  be  ordered  in  some  single  one-dimensional 
sequence  and  denoted  by  <l/)^  >  where  a  =  1,2,3,U,5,  ...  is  a  serial  index  whose 
values  label  individual  variables.  Then  for  given  k  the  realizable  values  of 
u(k),  u(-J<),  w(k),  w(-k)  are  fixed  completely  by  the  values  of  eight  of  the  q^. 
The  a  are  analogous  to  the  canonical  coordinates  and  momenta  of  Hamiltonian 
systans,  as  we  shall  see. 

It  follows  from  the  linear  relations  between  the  q's  and  the  u(k),  w(k) 
that  (3.1),  (3.2)  may  be  written  in  the  combined  foiro 

(3.7)  fia^W"    |:^  V^P  S-  *  ^a'     ^^^pr  '  W^' 

where  the  coefficients  A  _  and  damping  factors  v  are  functions  of  their  indices 
(but  not  of  the  q's  or  of  time),  and  the  f's  are  related  to  P^.(k)F.(k)  as  the  q's 
are  to  u.(k),  w.(k).  (VJe  shall  frequently  denote  the  row  vectors  with  elements  c^ 
and  f  by  q,  f.  No  summation  convention  will  be  used  for  the  Greek  indices.)  It 
is  readily  verified  from  (3.1),  (3.2),  and  the  definition  of  the  q's,  that  each 
V  is  positive  and  proportional  to  either  v  or  v. 

It  follows  from  (3.5)  and  the  definition  of  the  q's  that 
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(3.8)  ^-^Z%^- 

a 

We  remarked  earlier  that  if  v,   v  and  F(k)  are  put  equal  to  zero  then  (3.1), (3. 2), 
(3.k)  yield  E  -  0,     Therefore,   setting  the  v^  and  f^  equal  to  zero  in  (3.7),  we  find 

(3.9)  EeXl^a^-Z       ^a^%%\'^»  (v^,   f^-0). 


Since  any  assignment  of  initial  values  to  q  represents  admissibie  velocity  and 
magnetic  fields,  the  A    must  be  such  that  (3,S 
the  q's  equal  to  zero  except  q  >q'\  ,q_  »  we  find 


magnetic  fields,  the  A         must  be  such  that  (3.9)  holds  for  all  q.  Setting  all 


(3.10)         A  ._    +  A-  ♦  A_    . 

[i\a-         \(T\i        cr\i\ 


0. 


If  the  triangle  interaction  of  the  modes  <l,»q5.><3o-    is  defined  by  the  terns 

^\\a-\  V  *   ^h.(r  \i.%- %.'  ^^(r\i\%.  \  ^^  *^®  respective  equations  of  motion  for 
q  ,q^,q^ ,  it  follows  immediately  from  (3.10)  that  this  elementary  interaction 
makes  a  vanishing  contribution  to  t  in  (3.9)   and  therefore  is  individually 

conservative. 

The  'rate  of  energy-transfer  from  wave  number  k'  to  wave  number  k  •  frequently  is 
defined  as  Q(k,k' )  -  i  i   [k«H*(k-k')5*(k')»u(k)  -  k.u(k-k' )s(k' )'U*(k)],  since, (v,F=0), 
(8/8t)  -s  u  (k)'^(k)|   -  ^. ,  Q(k,k«)    [cfref.  U] .     The  conservation  property  is 
here  expressed  by  Q(k,k')   +  Q(k',k)  »  0,     Actually,  the  representation  of  the 
complex  of  triangle  interactions  as  a  complex  of  conservative  pair  interactions 
is  not  well-defined.     There  is  a  complementary  teim  Q(k»k-k' )  for  each  term 
Q(k,k')  in  ^lui  ^(JIS'iL'^*     These  two  terms  each  depend  on  all  thr«e  amplitudes 
u(k),  u(k')>  H.(}S~)S' ^   ^^^  their  complex  conjugates),  and  it  is  rather  capricious 
to  call  Q(k,k')  the  rate  of  transfer  f rom  k,'   and  Q(k,k-k')  that  fr«njc-k'.     This 
is  exhibited  by  noting  that  Q(k,k')  «  Q(k,k')  +  K(k,k', |k-k' ] ),  where  K  is  any 
function  antisymmetric  in  each  pair  of  its  arguments,  may  validly  replace  Q(]5>X*^» 
since  ^(k,k')  +  ^(k',k)  -  0,  and  ^(k,k')  +  Q(k,k-k • )  -  Q(k,k»)  +  Q(k,k-k').     It 
may  be  noted,  parenthetically,   that  (3.10)  is  a  stronger  statement  than 
Q(k,k')  +  Q(k',k)  ■  0,  since  it  does  not  involve  averaging  over  the  vector- 
components  of  u(k). 
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Most  of  the  A  „     defined  by  (3.7)  vanish  beceuse  of  properties  of  (3.1)» 
opY 
(3,2),     Vfe  first  note,   from  (3.U)  and  our  assumption  u(0)  -  w(0)   =  0,  that  only 

pairs  of  Fburier  coefficients  such  that  k'   and  k-k '   are  neither  parallel  nor 

antiparallel  to  k  (nor,   hence,  to  each  other)   can  contribute  to  the  right  sides 

of  (3»l),   (3.2).     The  immediate  consequence  is  that 

(3.11)  A^^    ■  0,         unless  o,p,Y  are  all  different. 

apY 

We  may  also  see  from  (3.1), (3.2)  that  the  trisngle  interactions  connect  q's 

representing  real  and  imaginary  parts  of  Fburier  coefficients  in  only  certain 

ways.     If  an  even  number  of  the  three  indices  a,p,Y  refer  to  imaginary  part  q's 

then  A  „   ,Ao  ^ ,A     „  all  vanish.     The  coefficients  describing  the  non-vanishing 
apf     py^     Yttp 

triangle  interactions  obey  the  relations 

^^•^^^  ^cx  6  Y.   "  ^a  p,Y     "  "'a.p  Y     "  ^a^p.Y, 

r^r'i  r^i'r  i^r'r  i^^i'i 

Here  a   and  a   are,  respectively,  the  real  and  imaginary  parts  of  a  single 

independent  vector  component  of  some  given  Fburier  coefficient,  with  a  similar 

interpretation  for  the  other  indices. 

Now  let  us  arrange  all  the  allowed  wave  vectors  in  some  one-dimensional 

sequence  jC/^,\,  k/p\,  ...  ,  Since  the  sianmations  on  the  right  sides  of  (3.1), (3. 2) 

extend  over  the  single  vector  k ' ,  the  number  of  A  „  which  have  some  single  fixed 

index  a,  p  or  Yj  and  which  connect  modes  with  wave  vectors  below  a  member  Jc/j,>  of 

the  sequence,  is  proportional  to  N  as  N  -►  oo.  The   number  of  nonvanlshing  A  ^^ 

apY 

with  two  fixed  indices  (implying  fixed  values  of  some  pair,  and  hence  all  three, 
of  k,  k',  and  k-k')  is  a  small  number  independent  of  ^   in  the  limit. 

U,  Formulation  of  the  statistical  problem 

U.l  Gibbs  ensembles.   The  statistical  mechanics  of  the  dynamical  system  obey- 
ing (3.7)  may  be  formulated  in  terms  of  Gibbs  ensembles  in  close  analogy  to  that  of 
Hamiltonian  systems.  The  application  of  Gibbs  statistical  mechanics  to  turbulence 
has  been  discussed  by  Burgers L°-l,  Onsager'-'-',  Hopf  •-  J,  and  other  authors.  Vfe 
shall  present  briefly  some  results  of  this  approach  and  then  point  out  its  inadequacy 
for  the  fun  statistical  description  of  stationary  turbulence. 

The  system  described  by  (3.7)  has  an  infinite  number  of  degrees  of 
freedom,  and  it  is  not  possible  to  speak  meaningfully  of  its  equilibrium 
statistical  mechanics  (v^  -  0,  f^^^  »  O).  We  avoid  this  difficulty  by  truncating 
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the  system,  removing  from  the  equations  of  motion  all  a     (and  f  )  such  that  a  >  N, 
where  N  is  an  arbitrarily  large  number.  With  a  suitable  ordering  of  the  q's,  this 
may  be  regarded  as  the  introduction  of  an  infinite  ohmic  or  viscous  damping  for 
all  modes  vd.th  wave  number  higher  than  some  limit.  The  necessity  of  this  artifice 
is  related  to  the  fact  that  the  Navier-Stokes  equation  for  sero  viscosity  does  not 
have  nontrivial  families  of  stationary  solutions,  but  leads  to  the  formation  of 
surfaces  of  discontinuity.  If  N  is  taken  large  enough,  the  truncated  system 
provides  a  valid  representation  for  the  viscous  case,  since  we  may  assume  that  the 
excitation  of  sufficiently  high  wave  numbers  is  negligible.  In  the  rest  of  this 
paper  we  shall  deal  with  the  truncated  system. 
According  to  (3.7),(3»ll)  we  have 

a 

This  is  equivalent  to  the  Liouville  theorem  obeyed  by  conservative  Hamiltonian 
systems,  and  it  shows  that  the  motion  conserves  measure  in  a  phase  space  with 
the  a     as  Cartesian  coordinates,  so  that  for  X  °  "^  '^^   time-invariant  ensembles 
defined  by  density  p(q)  are  those  which  satisfy  the  equation 


?^it\5    '"^'^^^■°' 


By  (3.8),(3»9)j  P  "  p(E)  is  an  obvious  solution.  More  generally,  the  solutions 
are  arbitrary  functions  of  functions  T^^Cq)  such  that  the  relations  ^(^(q)  "  const, 
are  the  integrals  of  the  simultaneous  system 


According  to  the  theory  of  such  systems,  we  may  expect  that  there  exist  N-1 
independent  y...     This  implies  that  there  are  N-1  constants  of  motion,  including 
the  energy,  associated  with  the  conservative  equations  of  motion,  with  a  corres- 
ponding restriction  on  the  ergodicity  of  the  truncated  system. 

For  f  =  0  the  density  p(q),  prescribed  at  some  instant,  determines  the 
ccmiplete  structure  of  the  distribution  of  solutions.  The  density  at  later  times 
may  be  determined  by  Liouville 's  equation,  and  averages  of  any  function  of 
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simultaneous  values  of  the  q's  may  be  evaluated  directly  by  integration  over  p. 

Averages  of  functions  of  nonsiitiultaneous  values,    such  as  <  (L(t)qo(t')  >,   also  are 

fixed  by  p,  but  their  evaluation  is  not  so  direct.     For  the  example  given, 

qo(t')  could  be  expressed  as  a  Taylor  series  about  time  t,   the  various  derivatives 
P 
in  the  series  expressed  in  terms  of  q(t)  by  repeated  use  of  the  equations  of  motion 

(and  their  time-derivatives),  and  the  result  integrated  over  p. 

Since  nonsimultaneous  averages  are  of  great  interest  in  turbulence,   the 
lack  of  symmetry  of  the  Gibbs  method  in  this  respect  is  a  disadvantage.     Moreover, 
when  external  forces  act,   the  Gibbs  method  no  longer  determines  nonsimultaneous 
averages.     To  see  this,  we  note  that  when  ^  /  0,   (3,7)  and  its  derivatives  still 
may  be  used  to  ejq^ress  higher  derivatives  of  q  in  terms  of  lower,  but  this  proce- 
dure introduces  time-derivatives  of  f.     In  this  case,   the  evaluation  of 
<  ^liC't)  <3q(*')  ^  ^y  ^^  Taylor  series  method  requires  knowledge  not  only  of  the 
distribution  of  simultaneous  values  of  q,  but  also  of  the  full  simultaneous  joint- 
distribution  of  these  values  with  those  of  f  and  all  time  derivatives  of  f .     This 
points  up  the  desirability  of  employing  a  more  complete  statistical  description 
than  is  provided  by  the  Gibbs  ensemble,  since  in  stationary  turbulence  external 
forces  are  a  necessity,  as  we  have  formulated  the  problem, 

U,2     Time-function  distributions.     The  full  solution  of  the  statistical 
problem  evidently  requires  the  determination  of  a  physically  appropriate  distribu- 
tion functional  T[q(    ),f(    )]   of  vectors  q(    ),   f(    )  which  obey  (3,7).*     Since  T 
describes  the  entire  history  of  the  systems  comprising  the  distribution,  it  does 
not  obey  an  equation  of  motion  like  the  Liouville  equation  for  the  time-chsnge  of 
the  Gibbs  density.     However,   a  functional  equation  for  I  can  be  constructed  which 
expresses  necessary  and  sufficient  conditions  th?t  the  members  of  the  distribu- 
tion obey  (3.7).     If  we  make  the  definition 

(U.2)         Ljt)  =  4^(t)  *  v^q^(t)  -    Z.  \pY  qp(t)  V*^  -  ^a(*)' 
this  equation  may  be  stated  as 

(U.3)  <Z  ^a^*^  ^a^*'^  ^T  '  °» 

a 

where  <  >«  indicates  an  average  over  Y[q(  ),f(  ) J .  Ihe  necessity  of  (ii,3)  is 

evident >  and  it  is  sufficient  if  imposed  for  t»t'.  This  follows  from  the  fact 

that  the  left  side  is  a  sum  with  positive  measure  over  a  svm   of  squares  of  real 
-  , 

At  this  stage  we  do  not  regard  the  distribution  of  external  forces  as  necessarily 
prescribed  a  priori.  The  restrictions  on  Y  appropriate  to  a  particular  physical 
situation  may  not  best  be  expressible  as  conditions  on  the  f's  alone. 
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quantities,   and  can  vanish  only  if  each  of  the  squares  vanishes,  except  for  a  set 
of  zero  measure.     The  possible  inclusion  in  the  distribution  of  a  zero-measure  set 
of  nonsolutions  is  of  no  physical  significance,  for  they  cannot  affect  finite-order 
moments.     The  integral  of  (U.3)   (with  t  =  t')  over  an  arbitrary  time  interval  is  the 
condition  that  the  members  of  the  distribution  obey  the  equations  of  motion  in  a 
mean  square  sense  within  the  integration  interval. 

Since  L  (t)   contains  only  linear  and  bilinear  terms  in  the  variables,   it 
is  clear  that  (U.3)   involves  only  second-,  third-,   and  fourth-order  moments  of  the 
distribution  of  q  and  f .     Therefore  we  need  evaluate  only  these  moments  to  determlns 
whether  a  given  distribution  satisfies  the  statistical  equations  of  all  orders 
obtained  by  multiplying  (3.7)  with  various  functions  of  q  and  f  and  averaging. 
This  suggests  a  high  degree  of  indetemiination  of  T  by  the  equations  of  motion, 
and,   as  indicated  before,  we  should  expect  as  a  minimum  requirement  for  determina- 
tion of  the  distribution  something  equivalent  to  complete  specification  of  the 
distribution  of  f ,  together  with  additional  conditions,   such  as  stationarity,  or 
specification  of  the  distribution  of  initial  values  of  q. 

It  should  be  noted,  however,  that  not  every  assignment  of  values  to 
second-,  third-,  and  fourth-order  moments  is  realizable.      The  requirement  that  T 
have  everywhere  non-negative  measure  in  the  function  space  of  q  and  ^  implies 
that  the  distribution  average  of  any  everywhere  non-negative  functional  of  these 
vectors  must  itself  be  non-negative'-     .     These  conditions  severely  limit  the 
realizable  moment  values. 

The  functional  Y  is  difficult  to  express  explicitly  because  it  depends 
on  a  nondenumerably  infinite  set  of  variables,  q(t)   and  f(t)  for  all  t.     We  may 
avoid  this  difficulty  by  confining  attention  to  an  arbitrarily  large  finite  time 
interval.     Ihen  we  may  expand  the  time-functions  within  this  interval  in  Fourier 
series  with  coefficients  q(oc)),  ?(w)  and  work  with  the  corresponding  functional 
[q(    )>£(    )]•     The  condition  on  J  equivalent  to  (U.3)   is 


(0      a 

where  L  (co)  is  the  Fourier  series  transform  of  L  (t), 
a  0. 

We  may  note  that  a  formal  solution  of  (U.U)  is 

J  =     H    N(a,co)  6  (L^(co))  Jj;  , 
a, CO  -^ 
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where  Jt  is  a  functional  of  J'(    )  only,   and  the  N's  are  normalization  factors. 
One  might  attempt  to  work  from  this  expression  by  employing  representations  of 
the  6-fanctions,  but  we  shall  not  take  this  approach  in  the  present  paper. 

ii.3     Stationarity  and  homogeneity,     A  stationary  distribution  is  one  whose 
moments  are  all  invariant  under  time  displacement.     That  is  to  say,   its  moments 
depending  on  quantities  with  arguments  t, ,  tp,   ...,  t     actually  depend  only  on 
the  n  -  1  independent  differences  of  these  times.     We  shall  assume  that  distidbu- 
tion  averages  over  a  stationary  distribution  can  be  identified  with  long-time 
averages  for  a  repi^sentative  physical  system,  macroscopically  isolated  in  the 
dissipationless  case,  or  subjected  to  stationary  driving  forces  in  the  dissipative 
case.     We  may  reconcile  this  ergodic  principle  with  the  previously  noted  lack  of 
ergodicity  by  appealing  to  random  perturbations  which  we  presinne  destroy  the 
constants  of  motion. 

In  a  stationary  distribution  only  the  following  moments,   and  their 
■o-derivatives,   appear  in  (U.3)! 


"""hXPyCl)   -  <  q^(t+Ti)  q^(t+t)  q^Ct)  q^(t)  >, 
Fjjp(Ti)   -  <  fjjCt+Tl)   fp(t)  >, 
Gj^pdi)   -  <  f^j(t+t)  qp(t)  >, 

We  shall  denote  these  moment  matrices,  with  indices  suppressed,  by  R,  S,  T,  F,  G,  H, 
From  the  definitions,  and  stationarity,  it  follows  that  R  Ax,)   ■  ^(~^)»  with 
similar  relations  for  the  other  moments.  In  terms  of  these  moments  (ii.3)  has  the 
form        f 

(U.6) 


Z  sy»^-fi;^(0  *  vX(^)  -  2  Z  V{V^'^  -  Vprx^^^  *  W"-^] 
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where  Sy]n{_z(t)  =  (l/2)[z(-o)  +  z(-t)],  for  any  z(t).  The  corresponding  expression 
for  (U.Ii)  involves  the  Fourier  series  transforms  of  the  moments  in  (li.6)  (defined 
in  an  arbitrarily  large  time  interval). 

Homogeneous  distributions  are  those  whose  moments  alD  are  invariant  under 
spatial  displacement.  Under  the  cyclic  boundary  conditions  we  have  adopted, 
homogeneity  requires  that  Fourier  coefficients  associated  with  diffei*ent  wave  number 
pairs  k,  -k  and  k',  -k'  be  uncorrelated.  If  the  distribution  also  is  invariant 
under  spatial  inversion,  the  real  part  of  a  given  Fourier  coefficient  must  be 
uncorrelated  with  the  imaginary  part.  Under  the  stated  conditions,  it  is  possible 
to  choose  orthogonal  principal  axes,  in  the  plane  normal  to  k,  along  which  the 
vector- components  of  u(k)  are  uncorrelated,  and  the  same  m^  be  done  for  w(k)  and 
F(k).  (The  choice  of  axes  may  have  to  be  different  in  each  case.)  Our  definition 
of  q  and  f_   leaves  us  free  to  take  their  elements  as  components  along  anj'  orthogonal 
axes  in  the  normal  planes}  let  us  choose  them  along  the  principal  axes  just  defined. 
Now  we  note,  following  Chandrasekhar^  ,  that  (3.1)>  (3.2)  are  invariant  under 
change  of  sign  of  w,  so  that  there  should  exist  distributions  in  which  u(k)  and 
w(k)  are  uncorrelated.  With  our  choice  of  axes,  it  follows  that  with  all  the 
conditions  stated, 

(U.7)    R^p(t)  -  6^p  R^(t>),       V^  =  6^p  F^(^). 

In  the  rest  of  the  present  paper  we  shfH  be  concerned  with  distributions 
obeying  (U.7),  and  in  which  q  and  f  have  zero  means.  For  their  existence,  we 
appeal  to  the  invariance  of  (3.1), (3. 2)  under  spacetime  displacement,  spatial 
inversion,  and  change  of  sign  of  w(k).  The  essuiTiption  of  homogeneous  external 

driving  forces  is  physically  unrealistic,  as  discussed  further  in  Section  10. 
This  assumption  clearly  is  required,  however,  in  any  consistent  idealization 
of  turbulent  motion  as  a  rigorously  homogeneous  and  stationary  phenomenon, 

U«U  The  hypothesis  of  weak  dependence.  Let  us  suppose  that  the  forces  f  , 
which  maintain  a  state  of  stationary,  homogeneous  turbulence,  are  statistically 
independent,  or  nearly  so.  We  hypothesize  that  in  turbulence  so  maintained  the 
statistical  dependence  among  any  group  of  n  modes  is  weak  if  n  is  small  compared 
to  the  total  number  of  modes  significantly  excited.   More  precisely,  we  assume 


The  order  of  magnitude  of  the  total  number  of  significantly  excited  modes,  may 
be  estimated  very  cnidely  as  the  cube  of  the  ratio  of  macroscale  to  microscale. 
It  is  quite  large  even  for  moderate  Reynolds  numbers. 
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that  normalized  moments  formed  from  the  a     in  such  a  group  are  closely  approxi- 
mated by  the  values   they  would  have  in  an  independent  distribution,   if  the  latter 
values  do  not  vanish,  and  are  small  compared  to  unity  if  the  corresponding  values 
in  an  independent  distribution  do  vanish.     This  hypothesis  is  based  on  the  fact 
that  the  motion  of  each  a     is  determined  by  coupling  to  f     and  by  coupling,  in 
an  individual,  complicated  fashion,  to  many  other  excited  modes,  so  that  strong 
dependencies  among  any  small  group  of  modes  would  be  difficult  to  account  for. 
The  further  justificstion  of  our  hjrpothesis  must  lie  in  the  consistency  of  the 
results  we  shall  obtain  by  its  use,   and,  more  logically,  in  constructions,  or 
existence  proofs,   for  distributions  obeying  it,  as  will  be  discussed  in  Section  8, 

It  should  be  noted  carefully  that  we  may  not  assume  the  weak  cross- 
dependencies  to  be  zero;  there  are  very  many  of  them,   and  together  they  represent 
the  entire  effect  of  the  nonlinear  interaction.     In  particular,  we  may  not 

assume  that  the  curnulants  associated  with  fourth-order  cross-moments  of  the 

[lOl 
q's  vanish.     This  assumption,  which  has  been  used  by  Heisenberg'-    -'  , 

Chandrasekhar '-     ■' ,   and  others,   is  seriously  incompatible  with  the  equations 

of  motion^    ■' . 

We  shall  assiune  in  the  remainder  of  this  paper  that  distributions 

appropriate  to  the  description  of  stationary,  homogeneous  turbulence  exhibit 

weak  dependence  as  we  have  defined  it. 

5.  The  Perturbation  treatment 

5.1  The  perturbation  hypothesis.  Let  us  consider  a  modification  of  (3.7)  in 
which  the  triangle  interaction  of  q^,  q.  ,  q   is  removed,  the  coefficients  A  ,   , 
A,^  ,  A^^  (originally  not  all  zero)  being  placed  equal  to  zero,  with  the  other 
A's  unchanged.  Suppose  we  seek,  according  to  some  suitable  criterion  ,  the  most 
nearly  independent  distribution  of  the  o     and  f  consistent  with  the  modified 
equations  of  motion,  and  with  given  ^  (^)  (or  given  E  for  non-dissipative  equilibrium, 
v^,  f  «  O),  The  interaction  we  have  removed  represents  a  small  part  of  the  total 
interaction,  if  many  modes  are  excitedj  therefore,  the  most  nearly  independent 
distribution  should  not  differ  very  much  for  the  orirjinal  and  modified  problems. 
However,  we  note  that  the  moments  S  .   (t),  S.   (z) ,   S   .(t)  do  not  appear  at 
all  in  (U.6)  for  the  modified  problem.  Since  this  equation  is  a  sufficient  condi- 
tion for  T  to  satisfy  the  equations  of  motion,  our  immediate  surmise  might  be  that 


Such  a  criterion  will  be  presented  in  Section  8, 
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in  the  most  independent  distribution  possible  these  moments  must  certainly  be 
very  small,  or  zero,   since  their  values  cannot  affect  (U.6). 

Matters  are  not  this  simple,  because,  if  other  cross-moments  do  not 
vanish,   zero  values  of  S  •^^('^)f  S,        (t),  S       At,)  may  not  be  realizable  in 
a  distribution  with  non-negative  measure.     On  the  basis  of  the  weak  dependence 
hypothesis  it  seems  plausible  to  sssiune,   however,  that  the  effect  of  the  realizability 
conditions  in  a  suitably  defined  most  nearly  independent  distribution  may  be  to 
induce  values  for  S  ,      (t),  S^       (x,),  S       Az)  at  most  of  small  magnitude  compared 
to  their  values  in  the  original  problem.     We  therefore  are  led  to  the  hypothesis 
that  in  the  original  problem  S  ^^(■'^)>  ^Xrru^^^*  ^(T   X^^^  ^^  "^■^•'"  appro'^J^ated 
by  the  contribution  to  their  values  obtained  by  introducing  the  triangle  inter- 
action with  coefficients  A  .         \cru»  ^cr  llX  ^^  ^  perturbation  on  the  motion  due 
to  the  rest  of  the  interactions. 

The  argument  given  may  be  extended  to  K(-c),   and  those  elements  of  T{z) 
which  vanish  in  a  completely  independent  distribution  of  q.     We  notice  that 
H  „  (t)  appears  in  (li.6)  multiplied  by  A  £.    ;  therefore  we  hypothesize  that  it 
is  well  approximated  by  the  contribution  to  its  value  obtained  by  introducing 
as  simultaneous  perturbations  on  the  motion  the  triangle  interaction  with  coef- 
ficients Ag„   ,   Ap     ,  A         and  the  force   canponent  f^.      (If  f,  is  removed  from  (3.7), 
H^g^(t)=0).    Similarly,  we  notice  that  the  elements  of  T(^)   with  indices  |j.,\,P,y 
in  some  order  appear  in  (U.6)  multiplied  bilinearly  with  A's  involving  these 
indices.     If  ia.,X,p,Y  are  equal  in  pairs,  these  elements  of  T(-c)  are  well  approxi- 
mated by  their  values  in  an  independent  distribution,  according  to  the  weak  de- 
pendence hypothesis.     Following  the  argument  above,  we  hypothesize  that  for  n,X,p,Y 
not  equal  in  pairs,   an  element  T  ^^^,^(1;)  is  well  approximated  by  the  contribution 
to  its  value  obtained  by  introducing  in  pairs,   as  simultaneous  perturbations  on 
the  motion,  those  tiriangle  interactions  whose  coefficients  bilinearly  multiply 
T  ^^.   (t)  and  the  elements  obtained  from  it  by  index  permutation. 

VJe  complete  this  statement  of  the  perturbation  hypothesis  by  noting 
the  evident     fact  that    the  element  G     (i)  is  given  by  the  contribution  to  its 
value  obtained  by  introducing  f     as  a  perturbation  on  the  motion. 

The  perturbation  hypothesis  just  stated  has  a  fairly  clear  intuitive 
meaning.     Consider  the  triple  conflation  of  a  ,  q^,  q     ,  described,   in  part,  by 
S  «      (t),  S,        (li),  S       •x^'^)*     These  three  modes  interact  botli  directly,   as  a 
result  of  the  terms  2A^^^q^  q^  ,  ^\^^<i^(i^f  ^^(TllX  \  \  ^°  *^®^^  respective 
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equations of  motion,  and  indirectly,  through  chains  of  other  modes  connected  by 
the  equations  of  motion.  Our  hypothesis  implies  that  the  triple  correlation  is 
entirely,  or  almost  entirely,  the  result  of  the  direct  interaction,  the  small 
contributions  of  the  individual  chains  of  indirect  interaction  combining  to  form 
a  small  or  zero  resultant.  This  certainly  seems  plausible  if  we  picture  turbulence 
as  a  mixing  process  in  which  phase  information  is  degraded  and  destroyed  by  the 
motion  as  a  whole.  In  similar  fashion,  it  may  be  seen  in  the  case  of  H  and  T 
also  that  our  hypothesis  picks  out  the  contributions  of  those  triangle  interactions 
which  constitute  the  most  direct  paths  of  interaction  of  the  modes  involved. 

At  the  end  of  this  section  we  shall  give  a  posteriori  support  to  the 
consistency  of  the  perturbation  hypothesis  by  actually  comparing  direct  and 
indirect  contributions  to  the  moments.  A  thoroughgoing  proof  of  validity  would 
depend  on  the  construction  of  distribution  functionals  obeying  the  hj'pothesis, 
as  discussed  in  Section  8, 

The  moments  defined  in  (U,?)  all  depend  on  values  at  only  two  times. 
The  arguments  given  above  can  be  restated  for  more  general  moments,  and  we 
therefore  extend  the  perturbation  hypothesis  to  third  -  order  moments  depending 
on  three  times,  and  fourth-order  moments  dependin?^  on  four  times. 

5,2  Perturbation  expressions  for  skew  moments.  We  shall  call  skew  moments 
those  cross-moments  which  must  vanish  in  a  completely  independent  distribution 
of  the  a  and  f  .  In  order  to  apply  the  perturbation  hypothesis  to  the  evalua- 
tion of  the  skew  moments  in  (h.6)  we  first  consider  the  effect  on  the  motion  of 
q  produced  by  introducing  at  t  »  0  an  arbitrary  perturbation  ?(t)  on  the  right  side 
of  (3.7),  If  K   is  not  too  large,  and  t  is  not  too  large,  in  a  sense  to  be  dis- 
cussed later,  th^  response  can  be  written 

t  t  t 

(5.1)  6q^(t)  "21   rs:j^^(t,s)C^(s)ds  +  Z  f  C  ri^^^{t,3,T)K^{8)l^{T)d5dT   + 


•J 


'o  *  cf    o 


where 


)C  a(*iS>  «  0,  (t  <  3),     rj^    ^^(t,s,r)  =0  (t  <  s  or  r),  ...  , 


and  6q  (t)  is  the  deviation  of  q^(t)  from  its  unperturbed  value.     In  general, 
K  ^(t,s),  »  ,      (t,s,r),    ...   are  functionals  of  the  unperturt)ed  q^  and  ^  and,   in 
addition,  depend  explicitly  on  their  arguments.     For  a  linear  system,  only  the 
first  term  of  this  series  is  non-zero,   and  ^T,(t,s)  is  simply  the  iirpulse  response 
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matri-X  of  the  system.  In  the  nonlinear  case  at  hand,  "C   ,(t,s)  gives  the  response 
to  an  infinitesimal  impulse. 

Now,  in  order  to  evaluate  the  triple  moment  of  a  ,  q^,  q,^  by  the  perturba- 
tion hypothesis,  we  take  K^   -  SA^^^q^q^,  K^   =  ^^^1f%.>   ^  '  ^yn^%%>   *° 
represent  the  switching  on  of  the  triangle  interaction  of  these  modes.  To  first 
order  we  have 

(5.2)      <  c^(t-L)qp(t)q^(t)  >  --  <  6q^(t--c)qp(t)q^(t)  > 

+  <  q^(t--o)5qp(t)q^(t)  >  +  <  <^(t-t)qp(t)6q^(t)  >, 

where  6qL(t-x),  5qg(t),  6q,^(t)   are  given  by  the  first  term  in  series   (5.1), 

The  weak-dependence  hypothesis  implies  that  ^       is  at  most  weakly 

statistically  dependent  on  q^  and  q^.     A  further  consequence  of  the  hypothesis  is 

that  on  the  average,  ^  „  and  J^       are  very  small  compared  to  the  diagonal  element 

ap  ay 

J^     ,  since  J^„o(t,s),  to  talce  one  of  the  pair,  is  the  response  of  q^  to  a  perturta- 

QO  up  Tl 

tion  in  the  equation  of  motion  of  q^,  to  which  it  is  only  weakly  coupled.  It 
follows  that  the  first  term  on  the  right  of  (5,2)  is  well  approximated  by 


(5.3)    <  5(^(t--c)qp(t)q^(t)  > 

t--c 
^  2A^PY  j    <  i^j^(t--o,s)  >  <  qp(t)qY(t)qp(s)q^(s)  >  ds 


o 

t-t 
2A^p^  (    <^„„(t-^»s)  >  Rao(t-s)R_(t-s)d3, 


o 


with  corresponding  expressions  for  the  other  two  terms.  The  averages  are  over 
the  unperturbed  distribution,  and  we  have  used  (3.1l)j  in  addition  to  the  weak 
dependence  hypothesis,  to  obtain  the  third  member. 

Now  we  note  that  since  the  unperturbed  distribution  is  stationary, 
<  X.   ■v(t,s)  >  can  depend  only  on  the  difference  of  its  arguments,  so  that  we  may 
define  the  distribution-averaged  infinitesimal  impulse  response  matrix  by 

(5.U) 


g^^(t-s)  .  <  ^^^(t-s)  >, 


g^^(t-s)  »  0,   (t  <  s). 
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(We  shell  refer  to  g"  as  simply  the  impulse  response  matrix.)  We  may  now  rewrite 
(5.3)  as 

(5.5)    <  5q^(t--o)qp(t)q^(t)  >  -  2A^p^  j   g^^(-t-3)Rpp(-s)R^(-s)d8. 

-t 

We  should  note  that,  for  a  given  pertuitation  strength,  we  cannot  expect 
series  (5.1)  to  converge  for  arbitrarily  large  t.  In  general,  we  should  expect 
nonlinear  response  effects,  represented  by  the  second  and  succeeding  terms  of 
the  series,  to  become  of  importance  when  6q(t)  predicted  by  the  linear  term 
becomes  comparable  with  typical  unperturbed  values  of  q.  At  Reynolds  numbers 
characterizing  fully  developed  turbulence,  a  given  triangle  interaction  represents 
only  a  small  part  of  the  total  force  acting  on  a  given  modej  therefore  we  may 
expect  the  error  due  to  approximating  series  (5.1)  by  its  first  term  to  be 
significant  only  for  values  of  t  of  the  order  many  times  the  characteristic  period 
of  the  mode  whose  perturbation  by  this  interaction  we  ai^  considering. 

Returning  to  (5.5^^^  we  note  that  RpgC^i)*  ^^^^^  ™^st  fall  off  with  increase 
of  z   rapidly  enough  for  j  R|,g(t)dt,  [  R  (x)dt  to  be  finite,  and  we  may 


anti- 
cipate that  g  (t.)  goes  to  zero  at  z   '="oo.*-  It  then  is  clear,  providing  z   itself 
is  not  large  compared  to  characteristic  periods,  th^t  the  contribution  to  (5.5) 
from  large  argrments  -z-s   is  small,  and  that  for  large  t  we  may  replace  the  limit 
-t  by  -00  without  appreciably  altering  the  value  of  the  integrsl.  The  new 
expression  depends  only  on  the  difference  time  tj  it  therefore  corresponds  to  the 
new  stationary  state  which  eventually  is  reached  after  switching  on  the  perturba- 
tion. For  z   large  compared  to  characteristic  periods,  the  first-order  perturba- 
tion expression  may  involve  a  si.giificant  fractional  error,  but  for  such  z   the 
triple  correlation  is  very  small  in  any  event. 

Working  out  for  the  other  two  terms  in  (5.2)  the  expressions  corres- 
ponding to  {$S)»   setting  t  =  oo,  and  using  R  (s)  =  ^uyS~^^ *   ^®  obtain  the 
complete  first-order  perturbation  result: 

This  is  reasonably  clear,  aside  from  other  arguments,  if  the  system  is  dissipative, 
For  the  nondissipative  case,  see  Section  6, 
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(5.6)         Sp^(t)  =  2A^p^  j    g^(s-t)Rpp(s)R^(8)ds 


.00 


*  2Ap^         ^aa(^-^)epp(^^V^"^^" 


o 

.00 


{■■ 


+  2A_-         R  „(s-t)R3o(s)B_(s)ds. 


o 

Next  we  may  evaluate  ^     (^)  by  introducing  f     as  a  perturbation  in  the 
equation  of  motion  of  a  .     At  first  sight  it  appears  thst  first-order  perturbation 
theory  should  be     inapplicable,   since,  at  least  for  modes  in  the  energy-containing 
range,  f     can  hardly  be  regarded  as  a  small  perturbation.     Actually,  we  may  use 
the  first-order  theory  if  f     has  a  normal,  or  nearly  normal,  univariate  distribu- 
tion.    In  this  case,  we  may  consider  f     the  sum  of  a  large  number  of  small  independent 
increments,   and  evaluate  the  covariance  of  a     with  each  of  these  increments  by  the 
fir^t-order  theory.     Writing  f  _   "  ^  a  f     ,  we  obtain  in  direct  analogy  to  the 
preceding  analysis 

CO 

where  F  nTi(s)  =  <  ^an^**s)f  (t)  >.  Summing  over  all  the  f^^,  we  have 

00 

(5.8)  G^Jx)  .  j"    §^(s-r)F^^(s)ds, 

n 

since  F (z)  =  T"  a  F         (z),  by  virtue  of  the  independence  of  the  f.     This  is 

aa  ^—     n  aann     ' '     "  ^^  an 

precisely  the  resvilt  we  would  have  obtained  by  considering  f     itself  a  small 
perturbation. 

The  evaluation  of  H  q.  (t)  and  T  ,<,   (t)  according  to  the  perturbation 
hypothesis  entails  the  use  of  second-order  perturbation  theory,   since  contribu- 
tions bilinear  in  pairs  of  couplings  are  required.     The  second-order  perturbation 

corresponding  to  Vi^^iz)  is 
apY 

(5.9)  62  ^  fa(t*^)qft(t)qy(t)  >  -  <  fjj(t+x)5^qp(t)6^q^(t)  > 

+  <  fjj(t-Hi)62q  (t)q^(t)  >  +  <  fi^{\.^z)c^{\.)b^{\)  >. 


-  19  - 


The  first-oixier  perturbations  ^^Qo   and  6-,<l^  may  be  evaluated  as  before.      The 

second-order  perturbations  ^2%*  ^2%  °°"sist  of  contributions  arising  frcmi  the 

second  term  of  series  (5.1)  and  of  contributions  obtained  by  substituting  for  K 

in  the  first  term  of  (5.1)  its  first-order  perturbation  6.5.     In  the  case  at  hand 

J.—' 

the  first  order  perturbations  indicated  by  the  perturbation  hypothesis  are 

etc.  I 

It  is  clear  that  the  leading  contributions  to  (5.9)  from  the  second  term  in 

series  (5.1)  involve  distribution  averages  of  products  of  odd  numbers  of  q's.  On 

the  basis  of  weak  dependence,  this  contribution  may  be  neglected  in  comparison  to 

that  arising  from  double  application  of  the  first  order  perturbation  formulas, 

which  gives  expressions  with  an  even  number  of  q  factors.  By  fiypothesis,  H  r^(t) 

is  well  approximated  by  those  contributions  arising  from  perturbations  of  the  q's 

which  involve  both  f  and  one  of  the  coefficients  A„o.,  Aq„,  A^q.  This  excludes 

a  'ipY   pY^  y^P 

the  first  term  on  the  right  side  of  (5.9).  Then,  evaluating  5-qp(t)  and  f'_q^^(t) 
as  was  indicated,  retaining  only  teims  not  excluded  by  the  perturbation  hypothesis 
and  the  weak  dependence  hypothesis,  we  find 

00    00 
(5.10)    H^^(x)  -  2  I  ds  I  dr  g^(r-8)[Ap^gpp(3)R^(s) 


F„-(r+-c). 


*   A^pg^(s)Rpp(s)|  -^ 


It  is  reassuring  as  to  consistency  that,  in  arriving  at  (5.10),  the  terms  excluded 

by  the  perturbation  hypothesis  a]*e  those  which  are  small  anyhow  as  a  consequence 

of  weak  dependence. 

According  to  the  weak  dependence  hypothesis,  the  non-skew  cross-moments 

Tg  ^(-c)  appearing  in  (U.6)  are  well  approximated  by  expressions,  in  terms  of  R, 

corresponding  to  an  independent  distribution.  The  evaluation  of  the  skew  Tq_^  .(t) 

proceeds  in  a  fashion  very  similar  to  that  for  H  „  (t).  Thei^  are  many  second- 

order  terms  involved,  however,  and  the  result  is  complicated.  It  may  be  seen 

from  the  selection  rules  on  the  A's  that  for  all  of  the  skew  '^OY^,^(''i)  except  a 

fraction  proportional  to  l/H,   where  N  is  the  total  nvmber  of  modes  in  the  truncated 

system,  it  is  the  case  that  if  T^  .(x)  appears  in  ([1.6)  then  there  are  no  non- 

PY|x\ 

vanishing  A's  with  one  index  p  or  y  and  another  index  simultaneously  |i.  or  X.  With 
these  conditions  (which  result  in  a  simplification). 
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a  o  ^ 

+  3  similar  terras  obtained  by  permuting  3  with  y  and/or  n  with  X, 

+  a  similar  term  obtained  by  permuting  p  with  y> 
rOO         >.oo 


C  "C  s 

+  a  similar  term  obtained  by  permuting  p,  with  \, 


(PjY»M.A  all  different). 
The  moments  given  by  (5.11)   are  not  zero,  in  general.     Thus,  second  and  fourth 
moments  are  not  related  as  in  a  normal  distribution.     The  expressions  for  the 
remaining  skew  Tg.   .  (t)  are  similar  to  (5.11). 

We  shell  now  present,  in  skeleton  form,  an  argument  in  support  of  the 
£  posteriori  plausibility  of  our  perturbation  h3rpothesis  for  the  calculation  of  the 
moments  S-      {z) ,     Similar  arguments  exist  for  the  other  moments.     An  evident  check 
on  the  consistency  of  the  hypothesis  is  the  comparison  with  (5.6)  of  the  perturbation- 
theory  contributions  to  Sc.^(^)  from  successively  longer  chains  of  indirect  inter- 
action  of  a  ,  q_,  a  ,  through  triangle  interactions  with  other  modes.     We  shall 
indicate  the  essential  results  of  such  a  program.     First,   it  is  not  hard  to  verify 
that  contributions  from  chains  of  even  numbers  of  triangle  interactions  are  negli- 
gible.    Crudely  speaking,   this  is  because  such  chains  lead  to  distribution  averages 
over  products  of  odd  numbers  of  the  q's.     The  leading  terms  in  the  contributions 
fixjm  chains  of  odd  numbers  of  tin. angle  interactions  cancel  in  pairs  as  a  result  of 
symmetry  properties  of  the  A's.     For  example,   it  is  not  hard  to  verify  that  typical 
leading  third-order  contributions  (involving  a  chain  of  three  indirect  interactions) 
are  of  the  form    ^      A^  ^A,       A         Z(p,Y,a,n,X, cr  )  where  Z  is  a  functional  of  the 

diagonal  elements,  with  the  indices  indicated,  of  the  matrices  g(    )   and  R(    ),     For 
every  terra  in  this  summation  there  is  a  conjugate  term  A     i>|A,  ,      ,A    .   ,  Z(p,Y>a,n',X',<T') 
in  which  one  of  the  pair  a  ,  a   ,   is  the  imaginary  part  of  the  Fourier  coefficient  of 
which  the  other  is  the  real  part,  with  similar  meanings  for  the  other  primed  indices, 
A  homogeneous  distribution  is  necessarily  symnetric  with  respect  to  real  and  imaginary 

parts  in  the  sense  g     (t)  "  g„.Mt(^)»  ^uu^^^  "  Vv*^^^'   ®*°*     Therefore 
Z(p,Y,a,|i,X,o-)  -  Z(p,Y,a,|ji«,X',cr-').     However,  it  readily  may  be  verified  from  (3.12) 
that  for  any  choice  of  p,Y,o,^,X,(r  we  have  A^^^A^^^A^^  -  "V'^'^'YO-^a" V'a* 
Thus  the  terms  in  the  summation  cancel  in  pairs.     The  same  result  may  be  verified  for 
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the  other  leading  thiivl-order  contributions  to  Sg^Ct)  in  which  the  indices  PjYi^i 
appear  differently.     The  order  of  magnitude  of  the  individual  third-order  terms 
is  such  that,   if  they  all  had  the  same  sign,   the  total  third-order  contribution 
would  be  the  same  order  of  magnitude  as  the  direct  contribution.     The  effect  of 
the  cancellation  is  to  insure  that  the  total  third-order  indirect  contribution 
is  negligible  compared  to  the  direct  contribution  (5.6),  if,   indeed,  it  does  not 
vanish  exactly.     The  same  is  true  of  the  other  highei^-order  contributions  to 
SQ^(t),  and  of  the  indirect  contributions  to   the  other  skew  moments,  to  the 
extent  that  they  have  been  investigated. 

There  exists,   therefore,  the  possibility  that  the  perturbation  hypothesis 
is  exactly  true  for  homogeneous  distributions.     The  more  rigorous  development 
of  the  argument  just  presented  is  most  appropriately  carried  out  in  terms  of 
explicit  distribution  functionals  of  the  type  discussed  in  Section  8, 

6,     Covariance  equations  for  non-dissipative  equilibrium 

When  V   ,f     "0,  and  the  simultaneous  distribution  of  the  q_   is  Maxwell- 
o     a  '  ^  ^ 

Boltzmann  (M-B),  a  very  simple  exact  expression  may  be  obtained  for  g(x;),  using 
only  (3.8)   and  (U.l)  and  the  fact  that  the  equations  of  motion  are  of  first 
order. 

Consider  two  conservative  systems  with  variables  a     and  p    which  obey  the 
equations  of  motion 

(6.1)  q^   -  C^(q,t),  p^  -  Pa^P'*^> 

have  energies 

a  a 

and  obey  the  Liouville  theorems 


a 

Consider  M-B  distributions  of  each  of  these  systems  with  identical  temperatures. 

The  joint  distribution  of  simultaneous  amplitudes  then  is  a  function  only  of  the 

total  energy  E^  +  E  ,  and  is  time- invariant, 
q   P 
Let  us  introduce  a  conservative  coupling  by  altering  the  equations  of 
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motion  to  the  form 

(6.U)     ^  -  <^  *  e  2:  A  (t)pp,   P<,  -  P,  -  e  Z  ^p^(t)q  , 

p  P 

whei^  A  rt(t)  is  an  arbitrary  matrix  function  of  time  which  vanishes  for  t  <  0, 
and  6  is  an  infinitesimal  parameter.  Equation  (6.3)  is  unaffected,  and  the  joint 
simultaneous  distribution  is  invariant  under  the  perturbation.  It  follows  that 

(6.5)     6  <  q^(t)p^(t)  >  -  <  6q^(t)p^(t)  >  +  <  q^(t)6p^(t)  >  -  0, 

for  all  \i   and  \,  Introducing  distribution-averaged  infinitesimal  impulse  response 
matrices,  as  in  the  last  section,  and  noting  the  rigorous  independence  of  the  full 
all-time  unperturbed  distributions  of  the  q-systems  and  p-systems,  we  may  write 
(6.5)  as 


(6.6)    21 


o 

where  h('c)  and  U(l)  are,  respectively,  the  averaged  infinitesimal  impulse  response 
matrix  and  covariance  matrix  of  the  p's.  Equation  (6,6)  holds  for  arbitrary  t, 
arbitrary  A  „(s),  arbitrary  forms  of  R  _(t)  and  II^Q(^),  and  for  all  \i   and  \. 
This  is  possible  only  if 

It  follows  from  the  definition  of  g^,(t)>  and  the  form  of  (6.1),  that 
g  (-KD)  «  6  .  The  solution  of  (6,7)  is  then 

(6.8)  g^(^)  =  R^(t)/R^(0),  h^p(t)  »  \p(^)^xx(0),    (t  >  0), 

where  we  note  that  in  the  M-B  distribution  we  have  taken,  R  (O)  ■  6  „R  (O), 

\6^^^  ""  \b^XX^^^*  ^^^   ^  (0)  =  \x(^^»  ^°^  ^^   index  values.  In  matrix  notation 
our  final  result  for  the  q  system  is 

(6.9)  g(Ti)  =  R"^(0)R(t),    (t>0), 
g(t)  =0,  (l  <  0). 
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Equation  (6»9)  is  an  exact  result  for  an  M-B  distribution.  This  distribu- 
tion seems  the  most  appropriate  physically,  but  it  is  of  interest  to  note  that  for 
systems  with  very  many  degrees  of  freedom,  and  weak  dependence,  (6.9)  should  be  an 
excellent  approximation  in  a  simultaneous  distribution  which  is  any  reasonable 
function  of  E  .  To  see  this,  we  may  consider  introducing  an  arbitrary  infinitesimal 
conservative  coupling  not  with  a  second  system,  but  among  the  q's  themselves.  We 
then  obtain  (6.6)  with  h  and  U  replaced  by  f  and  R.  We  then  may  argue  that  A^d(s) 
may  be  chosen  to  give  arbitrary  coupling  between  modes  triiose  covariances  are 
essentially  unrelated,  and  we  are  thereby  led  to  the  equivalent  of  (6.7)j  and  to 
(6,9).  We  m,ay  note,  also,  that  for  a  linear  system  the  impulse  response  must  be 
independent  of  distribution  so  that  (6.9)  gives  g(t)  for  an  individual  system, 
if  R(t)  refers  to  any  distribution  displaying  equipartition.   (The  system  obey- 
ing (3.7)  is  not  linear,  of  course.) 

If  we  multiply  (3.7)  by  non-simultaneous  q^  and  average  over  i;  we  obtain 


(6.10)  R^(^)  -Z  A^p^Sp^(t),      (v^,  f^-  0). 

Expressing  equipartition  at  eouilibrium  by  R  (O)  -  R,  we  find,  using  (3.10 )j (5.6)  and 
(6.9),  ^  ^ 

(6.11)  Sp^(.)  -  -  ^  I  R^Jx-s)Rpp(s)R^(s)ds. 


We  see  that  S-,  (t)  is  an  odd  function  of  x.  (Sq (O)  must  vanish  since  the  Gibbs 

ensemble  is  a  function  only  of  E  '  -^  ^  q^  ,)     We  also  note  that  S-^(t)  goes  to 


«* 


lero  as  t  -»►  00,  provided 'R(tj  falls  off  properly  for  large  argument.   By  (6.10) 


*For  a  linear  system,  (6.9)  is  anticipated  by  the  work  of  H.B.  Callen  and  his  co- 
workers, who  use,  howefver,  a  different  appro ach L-'-^-' *<- ^J  .  The  relation  between 
the  correlation  and  response  matrices  associated  with  small  fluctuations  from 
equilibrium  of  macroscopic  thermodynamic  variables  is  examined  by  Callen  and  Greene 

An  interesting  estimate  on  the  order  of  magnitude  of  the  triple  moment  given  by 
(6.11)  may  readily  be  obtained  for  an  idealized  system.  Suppose  the  coupling  is 
such  that  all  the  modes  have  essentially  the  same  characteristic  period,  assume 
all  the  nonvanishing  A's  to  be  of  the  same  order  of  magnitude  and  let  the  average 
number  of  nonvanishing  A's  with  one  fixed  index  be  M.  Then  this  period  may  be 
estimated  by  squaring  (3.7)  and  averaging  [see  [l2l  for  the  methodj ,  and  conse- 
quently it  may  be  found  that  the  maximum  values  attained  by  the  normalised  triple 
moments  are  of  order  }C  '    ,     Since,  with  the  properties  given  in  (3.11),  seq.,  M 
is  the  order  of  the  total  number  of  excited  modes,  this  estimate  supports  the 
consistency  of  the  first-order  perturbation  treatment. 


Li5j 


-  2U  - 

and  (6.11), 

(A   )^  ^ 
(6.12)     k^M  *  2     Z  -^  (\a^^-^^   V^"^  V^^^  ^"  "  °- 


Noting  that  g^-C+O)  "  1>  we  may  rewrite  (6.12)  in  the  'universal'  forni 

X   y 
(6.13)    g^Jx)  =  1  -  2  X:  (A^pY^^  I  dy  f  dz  g^^(y-z)  g3p(«)  g^(z),  (x  >  0), 

1/2 

where  x  ■»  -cR  '  .  Assuming  the  homogeneity  conditions  expressed  by  (U.7)>  we  havB 

now  determined  the  complete  covariance  matrix,  for  (6.13)  niay  be  solved  by  step-by- 
step  numerical  intef:ration,  starting  at  x  =  0. 

We  may  verify  easily  that  if  the  moments  we  have  obtained  are   realizable 
they  satisfy  the  fundamental  statistical  equation  (li.6).  As  discussed  in  Section  h, 
(U.6)  is  satisfied  for  all  t  by  any  realizable  distribution  which  satisfies  it  for 
t  »  0.  For  V  ,  f  ,  ^  «  0,  (li,6)  becomes 

(6.Ui)   T  -R  (0)  +  2  ^  A  „  S-   (0)  +  T     A  „  A  .  T  .^  (0)  »  0. 

Now  T  .  -  (0)  vanishes  in  a  simultaneous  distribution  depending  only  on  E  =  r  ^  a  , 
unless  the  indices  are  equal  in  pairs.  If  we  assert  weak  statistical  dependence 
of  the  modes,  we  have  for  the  latter  case, 

(6.15)    T  ,„  (0)  ^  (5  «5^   +5  6--,  ♦  6  >6o  )R^» 

(n,X,p,Y  not. all  equal), 
so  that,  taking  account  of  (3.11), 

,  Z    ^  A       A^     T         (0)   ^   2    Z  (AapY)'  H^. 

This  is  an  exact  result  if  the  Gibbs  ensemble  is  M-B.     By  differentiation  of  (6.11) 

For  a  linear  dynamic  system,   the  equations  analogous  to   (6.13)  become  a  ssrstem 
of  simultaneous  bilinear  algebraic  equations  under  Laplace  transformation. 
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(noting  R  (O)  =  O)  we  obtain 


aa 


From  (6.16)  and  the  derivative  of  (6.10)  for  t  «  0,  we  obtain  (6,lU).  It  may  be 
verified,  using  the  perturbation  expressions  for  T(t),  that  (U.6)  is  satisfied  by 
our  moments  for  t  /  0  also, 

7.  Covariance  equations  for  the  dissipative  stationary  state 

Stationary  turbulence,  with  viscosity  and  driving  forces  acting,  represents 
a  state  of  strong  departure  from  equipartition  and  detailed  balance  in  which  the 
energy  of  the  system  fluctuates.  In  this  case  the  method  used  in  the  last  section 
to  evaluate  g(t)  is  inapplicable.  We  can  deduce,  however,  a  set  of  integro- 
differential  equations  which  simultaneously  determine  the  g-^(t)  and  the  R(,jj(t). 

Multiplying  (3.7)  by  nonsimultaneous  f  ,  and  averaging  over  T,  we  obtain 

(7.1)  -G^(t>)  -  v^G^(.)  .  ^  Vr^aPY^^^  '  \a^^^' 

p,Y   •^' 

If  now  we  introduce  an  infinitesimal  perturbation  6f  ,   statistically  independent 
of  f  ,  and  with  auto-covariance  5„F     (l),   the  perturbations  in  G     (t>)  and  H  Q^(t) 
may  be  obtained  immediately  from  (5. 8), ($,10)  by  substituting  S^^oa  ^°^  ^ao*     ^^*^ 
some  manipulation,   the  use  of  partial  integration  on  ^oG     ("c)   (noting  gL^C"*^)  "  1 
and  taking  5„F     (t)   zero  for  z  =oo),  we  obtain  the  following  expressions  for  the 
perturbations  in  the  terms  of  (7.1) J 

(7.2)  t^^iz)  s  62F^(t), 

o 

*  Vp  V^^^   V^'^]  1^^-'^  Vaa^'"^^'^^^* 
where    n(x)  -  |(x+|x|  )/|x| .     Equation  (7.1)  can  hold  for  arbitrary  e^F^^  only  if 
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the  coefficients  of  ^n^aa.^^^^   ^"  ^^^   integrands  of  its  seversl  terms  total  to 
zero.  We  thereby  find 

z 
(7.3)    ^(.)  -  v^g^(.)  =  2^  A^p^  f  ^,(-'^)[^^^S^^(s)R^is) 

Multiplying  (3.7)  by  nonsimultsneous  o  ,  and  averaging  over  T,  we  obtain 

(7.M   R^(t)  *  v^R^Cx.)  =  Z  V  V^'^  *   °aa(^) 

Equations  (7.3)»  (7.1;)>  (5.6),  and  (5.8)  form  a  complete  simultaneous  system 
which,  with  the  boundary  conditions  &,(,(+0)  "  Ij  '^on^'^^  "  ^'  determine  both  the 
g^^(t)  and  the  R„„(l)  in  terms  of  the  external  force  covariance  F  (z) ,     To 

uu  uu  QO 

verify  the  consistency  of  these  results  with  those  of  the  last  section,  we  take 
v^  "  ^aa   ■  0  i"  (7.3)>  {l»h)j   assert  equipartition,  and  assert  (6.9).  Then, 
using  (3»10),  we  readily  find  that  (7.3)  and  (7.U)  become  identical  with  each 
other,  and  with  (6.12),  It  may  be  verified  that  if  (7.3),  (7.U)  are  satisfied, 
then  (li.6)  is  satisfied  by  the  perturbation  theory  moments. 

At  first  sight  it  might  be  thought  that  the  system  (7.3)>  (7.U)  does  not 
constitute  a  much  more  difficult  problem  than  (6.12),  aside  from  the  increase  in 
the  number  of  unknowns  by  a  factor  of  two.  Actually,  the  difference  in  complexity 
is  profound.  Equations  (7.3) »  (7.ii)  cannot  be  solved  by  straightforward  step-by- 
step  integration,  starting  with  gLQ(*^)»  ^oa^^^*  ^^   *^®  first  place  the  depend- 
ence of  R_„(0)  on  a  is  one  of  the  answers  we  seek,  and,  moreover,  the  right  sides 
of  (7.3) »  (7.M  even  for  small  z   depend  on  g   and  R  for  all  argument  values.  The 
increase  in  complexity  is  not  surprising  when  we  consider  that  the  equilibrium 
nondissipative  problem  could  have  been  solved  in  elementary  fashion  by  expanding 
R  (Ti)  as  a  Taylor  series  about  z  "  0,   expressing  the  derivatives  in  terms  of 
simultaneous  averages  by  repeated  use  of  the  equations  of  motion,  and  averaging 
over  the  time-invariant  simultaneous  distribution.  No  such  direct  procedure 
exists  for  the  dissipative  case. 

Although  in  the  dissipative  case  there  may  be  no  general  relation  between 
g  j^(t.)  and  Rjy,(t)  of  simplicity  comparable  to  (6.9),  we  may  conjecture  that  at 
least  qualitatively  correct  results  will  be  given  by  the  assumption 
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(7.5)        e^M^  Y^ e    °-  ,  (l>0). 

aa(0) 

We  note  that  with  v     =  0,    (7.5)   reduces  to  (6,9),  and  (since  R     (O)  »  O)  it 
implies  ^(+0)  +  Vjjg^j^(+0)  -  0,   as  required  by  (7.3).     For  v^  very  large,   (7.5) 
implies  e     (t)  '~    exp(-v  ^),    (t  >  0),  the  correct  resiilt  for  the  impulse  response 
of  a  system  controlled  by  linear  damping.     This  property  indicates  at  least  the 
qualitative  correctness  of  (7.5)  for  heavily  damped  modes  in  the  dissipation 
range  of  turbulence.     For  the  inertial  and  low  wave  number  ranges,  where  the 
factor  exp(-v  t)  may  be  ignored,   the  qualitative  correctness  of  (7.5)  seems 
plausible  on  the  intuitive  ground  that  the  time  it  takes  a  mode  to  forget  a 
perturbation  impulse  through  the  mixing  action  of  the  other  modes  with  which 
it  interacts  should  be  approximated  by  the  mode's  own  characteristic  period. 

8,     Variational  construction  for  the  distribution  functional 

We  have,   in  principle,  determined  moments  which  satisfy  the  necessary 
and  sufficient  statistical  condition  (U.6).     The  complete  perturbation  theory 
solution  of  the  statistical  problem  requires  the  construction  of  a  distribution 
functional  which  yields  these  second-,   third-,  and  fourth-order  moments  and  fixes 
the  higher  moments.     A  proof  of  the  realizability  of  such  a  functional  would,  of 
course,  be  the  most  satisfying  check  on  the  validity  of  our  method.     It  was 
indicated  earlier  that  the  distribution  functional  certainly  is  not  fully  deter- 
mined by  its  second-,   third-,   and  fourth-order  moments.     We  need  an  additional 
condition  to  pick  out  a  particular  distribution  functional  appropriate  to  the 
physical  situation.     We  shall  adopt,   for  this  purpose,   a  variational  procedure 
which  defines  and  determines  the  most  nearly  normal  distribution  consistent  with 
the  given  moments,   and  with  any  additional  physical  conditions  which  can  be 
expressed  in  terns  of  distribution  moments. 

The  variational  method  is  best  foiroulated  in  terms  of  the  transform 
functional  J.    We  therefore  shall  summarize  the  foregoing  perturbation  theory 


In  this  connection  we  should  note  that  in  a  physical  situation  the  driving 
forces  exerted  on  low  wave  number  modes  (whose  behavior  is  substantially 
influenced  by  direct  action  of  the  driving  forces)  will  exhibit  characteristic 
periods  comparable  to  those  induced  by  the  nonlinear  forces  exerted  by  the 
interaction  with  the  higher  wave  numbers.     This  is  true,   for  example,   in 
turbulence  supported  by  a  steady  shear  flow. 
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in  terras  of  frequency  transforms.     In  order  to  deal  with  a  denunerably  infinite 

set  of  variables  we  choose  an  arbitrarily  long  time  interval  T  and  define  the  series 

transform  ^(co)  of  any  time  function  z(t)  by 


(8.1)         z(t)  =  T-^  2  z(a3)e^*,  z(co)   -   (    z(t)e'*"^*dt. 


(The  transform  of  any  qviantity  will  be  represented  by  the  corresponding  quantity 
with  a  wavy  overbar.)  In  the  limit  T  ->•  oo,  th©  series  transform  defined  in 
this  way  goes  over,  without  further  normalization,  into  the  integral  transform 
defined  by  z(t)  «  (27i)~  j   "^(01)6   d»,  provided 'z(co)  exists  in  the  limit. 
The  transform  of  the  averaged  response  of  mode  qL  to  an  infinitesimal  perturbation 
h^(t)  is 

(8.2)  ^  q^(co)  »  i^^(co)  h^(co), 

where  gjj^('«i)  is  the  transform  of  g-,_(^)»  The  expressions  for  the  perturbation- 
theory  moment  values  become  simpler,  in  some  respects,  in  the  transformed  repre- 
sentation. For  example,  the  transform  of  (5.3)  is 

(8.3)  G^(co)  '\^i-o^)y^{o>). 

We  shall  choose  as  fundamental  arguments  for  J  the  real  and  imaginary 
parts  of  the  quantities  C3  (co),  f  (co)  for  co  >  0.     (Tine  corresponding  variables 
for  CO  <  0  then  are  fixed  by  reality  requirements.)     For  simplicity  in  notation 
let  us  denote  all  these  independent  real  and  imaginary  parts  by  the  elements  y 
of  a  row  vector  y,  so  that  J  «  J(y).     We  may  assume  from  the  physics  of  our 
problem  that  frequencies  above  some  arbitrarily  high  limit  will  be  negligibly 
excited.     We  may  therefore  cut  off  the  sequence  y    at  some  arbitrarily  high  n 
and  deal  with  a  finite  number  of  variables. 

Let  us  write 

(8.!i)        lij)  -  Ay)0(y), 

which  expresses  the  non-negativity  of  J,  and  otherwise  involves  no  loss  of 
generality.     We  define  the  most  normal  distribution  consistent  with  given 
fixed  moment  values  as  that  distribution  for  which  the  probability  amplitude  0 
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minimizes 

(3.5,   iJ-'/-V'^ 
J  (Tj^  dy 

under  the  integral  constraints  represented  by  the  fixed  moment  values,  and  the 
boundary  condition  0(oo)  =  0.  In  our  problem,  the  fixed  moments  are  all  the 
second,  third,  and  fourth  moments  of  the  y  ,  as  given  by  the  perturbation  theory. 
Therefore  the  0  so  defined  obeys  the  Schroedinger  equation 

(8.6)   V^We-Z  X^y^y^  -  Z  ,  '^nm^  ^nV^  "  ^,      ^r^is^JJl^sV  ' 
^         I        n,m  n,m,£    ^        n,in,Jl,s  J 


0, 


where  6  is  the  lowest  eigenvalue  for  which  the  eoruation  has  -   solution  vanishing 

at  infinity,  and  the  Lagrange  multipliers  X^^,  X^^^^  ,  X^^^^g  are  uetermined  a 

posteriori  by  the  specified  moment  values.  It  is  apparent  that  we  may  take  0  real 

without  loss  of  generality. 

It  is  readily  verified  that  if  the  specified  moments  had  consisted  only 

of  all  the  second  moments  (X  .=0,X  »  =0)»  the  J  defined  by  this  criterion 

rmt  rmj?s  -^ 

would  have  been  the  Joint  normal  distribution  yielding  these  moments.     (The 
Schroedinger  equation  for  this  case  describes  an  arbitrarily  oriented  multi- 
dimensional harmonic  oscillator. ) 

The  solution  of  (8.6)  in  a  space  of  an  indefinitely  large  number  of  dimen- 
sions would  be  impracticable  except  for  the  fact  that  the  skew  moments  represent- 
ing departure  of  the  distribution  from  multivariate  normal  are  small,  so  that  the 

individual  nondiagonal  potential  terms  \^i7^7^j^  ,  Vnfs^nV^^s  ^  ^^'^^  "^^ 
be  regarded  as  small  perturbations  (with  suitable  precautions  and  subject  to 
complications).     In  this  case,   it  may  be  seen  to  follow  from  the  hypothesis  of 
weak  dependence  that  0  satisfying  (8.6)   gives  a  nearly  independent-Gaussian 
distribution  for  any  few  modes  when  averaged  over  the  distribution  of  the  other 
modes.     Then  the  perturbation  solution  of  (8.6)  may  be  written 

'An  alternate  variational  criterion  is  the  minimization  of  Jdyfln(f/J|dy)/j|dy 

under  the  same  Integral  constraints.     This  provides  a  direct  analogy  to  the 
principle  of  equal  a  priori  probabilities  based  on  the  H-theorem,  but  it  must  be 
remembered  that  |  is  not  the  Gibbs  probability  density  and  the  appropriateness 
of  this  principle  is  not  self-evident  in  the  present  context.     Although  the  alter- 
nate criterion  leads  more  readily  to  explicit  expressions  for  f ,  that  pi-esented  in 
the  text  has  a  certain  intrinsic  interest  because  of  the  appearance  of  the  proba- 
bility an^Dlitude  in  a  classical  problem.     It  is  a  conseanpnce  of  the  weak  dependence 
property  that  both  criteria  lead  to  nearly  the  same  distribution;  the  difference 
does  not  seem  significant,  particularly  in  view  of  the  irreducible  arbitrariness 
in  assigning  statistics  to  the  external  driving  forces. 
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(8.7)   0(y)  =  JL    iP^il   (1  *  i  c    y^y^y^  )    H    ^]iL  ^^*   l^ratuyr^syt^u^^^Z^' 
n,m,i?  r,s,t,u 

where  'X  (y)  is  a  joint  normal  distribution  amplitude  and  the  N's  are  normalization 

factors.  The  c„^  »  ,   d.^.  .  and  the  parameters  in  7C(y)  must  be  chosen  to  give  the 
nm  X*   rs  uu  — 

perturbation  theory  moment  values.  A  sufficient  condition  for  the  realisability  of 
our  moment  values  is  the  existence  of  such  c's  and  d'g. 

We  shall  reserve  for  a  subsequent  paper  the  evaluation  of  the  parameters, 
the  discussion  of  consistency  questions,  and  the  method  of  use  of  the  distribution 
functional.  We  wish  to  point  out  here  that,  assuming  the  success  of  the  prcedure, 
the  accuracy  of  the  predicted  higher  moments  may  be  expected  to  become  poor  when 
the  order  of  the  moment  is  comparable  to  some  function  of  the  total  number  of 
degrees  of  freedom  a  which  are  excited.  This  is  because  our  underlying  assumption 
of  weak  statistical  dependence  cannot  be  asserted  for  such  high-order  moments. 

We  must  also  postpone  discussion  of  the  interesting  question  of  the  close 
relation  of  distributions  of  the  form  defined  by  (8,6)  to  the  M-B  distribution  of 
simultaneous  values,  in  the  non-dissipative  case,  and  to  the  condition  of  minimum 
entropy  production  in  the  dissipative  case. 


9.  Isotropic  turbulence 

9.1  Zero  magnetic  field.  We  shall  now  present,  in  terms  of  the  usual  spectral 
representation,  the  closed  equations  for  the  velocity  covariance  which  insult  for 
isotropic  nonmagnetic  turbulence,  under  the  simplifying  assumption  (7.5).  Without 
loss  of  generality,  we  may  take  the  external  force  field  divergenceless.  Then  we 
may  define  the  scalar  functions  U,  F,  and  G  by 

(9.1)      <  u^(K,t  ♦  Ti)u*  (k,t)  >  •  ^  Pij(k)  U(k,Ti), 
<  F^(k,t  +  t)F*  (k,t)  >  -  ^  P^j(k)  F(k,t), 

Rel*^  ^±^h^  *  ^^   ^h^^  ^}''    •?  ^i:J^iS^  G(k,t).  (Re(}      denotes    'real  part') 

These  scalars  all  are  real,   and  U  and  F  are  even  functions  of  t.       U(k,t)   and 


*. 


Rigorous  isotropy  requires  L-»-oo,  with  our  boundary  conditions.  We  reserve 
for  a  aore  detailed  paper  the  discussion  of  the  validity  of  the  passage  to 
conditions  (9.1),  and  the  associated  passage  to  the  Fourier  integral  repre- 
sentation used  in  (9.8). 
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F(k,i),  suitably  normalized,  are  the  spatial  Fourier  coefficients  of  the  defining 
scalers  of  the  velocity  and  driving  force  spacetime  autocovariances,  and  k  U(k,0) 
is  proportional  to  the  usual  one-dimensional  wave  number  spectrum  function. 
We  may  obtain  from  (3.1)   (with  w  =  O)   the  equation 

(9.2)  U(k,x)  +  vk^U(k,t)     -     S(k,T;)  +  G(k,i), 
where 

(9.3)  S(k,t)  =  Im  [k^  T  <  Uj(k«,t)  u^  (k-k',t)  u*(k,t-Ti)  >|  , 

and  Im  {  }  denotes  imaginary  part.  The  quantity  k  S(k,0)  is  proportional  to  the 
total  net  rate  of  energy  transfer  by  nonlinear  interaction  to  all  the  modes  of 
wave  nxunber  k.  Evaluation  of  S(k,i)  and  Q(k,-c)  in  terms  of  U(k,-c)  and  F(k,t) 
may  be  carried  out  directly  on  the  basis  of  the  preceding  analysis,  since  the 
u,(k,t)  are  linear  homogeneous  functions  of  the  a.     Adopting  (7.5)  we  have  the 
expression 

(9.U)   g(k,i)  -  ^^^  exp(-vk^t),     {z>0), 

for  the  diagonal  impulse  response  of  a  mode  of  wave  number  k.  According  to  the 
perturbation  and  weak  dependence  hypotheses,  then , 

(9.5)  G(k,x)   =  Re  i<  T^{^,   t  +  i)5u^*(ii,t)  >},  (t  ^   oo), 
where                                . 

(9.6)  3u*  (k,t)  -  j  g(k,t-3)  F*  (k,s)  ds. 

o 
Proceeding  as  before,  we  find 

00 

(9.7)  G(k,t)  -  j    ds  gjj^ygj  exp(-vk^3)  F(k,s+t) 

o 

The  evaluation  of  S(k,i;)  is  more  involved  because  of  complicated  angle- 
averages  which  occur.  The  result,  expressed  in  terms  of  a  Fourier  integral 

F(k,t)  and  u(k,t)  are  each  linear  combinations  of  several  independent  variables; 
(U.7)  ensures  that  no  consequent  cross  terms  enter  (9.7)  and  (9.8), 


representation,   is 
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00    00      k+k ' 


(9.8)     S(k,T;)  "  TTk  r  ds  j  k'dk-  (       k"  dk«  [  ^^S(k](;^"  ^  exp(-vk^8)U(k,8)D(k',s+T;) 

"'o    ''o       |k-k'| 

XU(k",s+x)  -  ^fe^]5)"^  exp(-vk'28)U(k,3-L)U(k«,s)U(k",  s)], 

where 

(^^^)    a(k,k',k")  -  (1  -  xyz  -  2y2z^), 

b(k,k',k")  =  2k"^k'(-2  +  xy  +  x^z  +  y^z  +  2z^  *  2xyz^), 

and  x,y,z  are  the  cosines  of  the  interior  angles  opposite  k,  k',  k",  respectively, 

in  a  triangle  with  these  quantities  as  sides.  The  quantities  a,  b  obey  the 

relations 

a(k,k',k"  )  =  a(k,k"  ,k')  >  0, 

^9«10)        k^(k,k',k"  )  -  k«^(k',k,k"), 

b(k,k',k"  )  +  b(k,k"  ,k»)  =  2a(k,k',kn  ), 

in  the  domain  for  which  they  are  defined.  The  integration  indicated  in  (9.8) 
actually  is  symmetric  between  k'  and  k",  since,  if  |k  -  k'l  <  k"  <  k  +  k',  then 
|k  -  k"|  <  k'  <  k  +  k".   As  a  consequence,  it  may  be  seen  that  (9.10)  expi^sses 
the  consistency  of  the  expression  for  S(k,0)  given  by  (9.8)  with  the  detailed 
conservation  property  (3.10).  For  many  purposes,  a  and  b  may  be  approximated  by 
very  simple  forms. 

Equations  (9. 2), (9. 7), (9. 8),  together  with  the  boundary  condition  U(k,oo)=0^ 
form  a  closed  system  determining  U(k,t)  in  terms  of  F(k,t),  the  driving  force 


See  footnote  following  equation  (9.1).  The  normalization  of  the  integral  trans- 

foiTOs  is  fixed  by  U(k,Ti)  «  J  U(x,t)e"^*-d-'x,  where  U(x,t)  -  <  u(y+x,  t+i)-u(y,t)  >, 

with  similar  definitions  for  G(k,t),  F(k,t)  when  interpreted  as  integral  transforms. 
This  implies  a  corresponding  renormalization  of  u(k,t),  F(k,t)  to  unit  volume. 

The  integrand  term  in  (9.8)  involving  a(  )  indicates  a  positive  contribution  to  S(k,0) 
from  all  triangle  interactions  among  modes  of  wave  numbers  k,k',k",  while  that  involv- 
ing b(  )  indicates  a  negative  contribution.  Equation  (9,10)  assures  that  the  contri- 
butions to  k  S(k,0),  k'  S(k',0),  k"^(k%0)  from  these  same  interactions  add  to  zero 
as  required  by  detailed  conservation.   The  factors  l/[J(k,0)  in  the  positive  term  of 
S(k,0)and  lA^(k',0)  in  the  negative  term  indicate  th?t  the  theory  gives  a  net 
energy-transfer  in  the  direction  tending  to  establish  equipartition  among  the 


•K^ 


Fourier  modes* 
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covariance.     The  detailed  discussion  of  these  equations  vrill  be  reserved  for  a 
subsequent  paper,  but  we  shall  indicate  very  briefly     here     the  consistency  with 
them  of  Kolmogorov '  s  similarity  considerations.    -'*'• -^  The  Kolmogorov  theory  of  very 
high  Reynolds  numbers  predicts  the  existence  of  an  inertial  range  of  wave  numbers 
in  which  dissipation  is  negligible  and  for  which 

(9.11)  U(k,C)  =  (kA^)'^/^  U(kpO), 

if  both  k  and  k,  are  in  this  range.  This  spectrum  form  follows  firmly  frorr 
dimensional  analysis  if  it  can  be  established  that  energy  transfer  in  the 
inertial  range  proceeds  by  a  cascade  process  through  a  chain  of  successively 
higher  wave  n\uiibers,  the  essential  requirement  being  that  no  appreciable  fraction 
of  the  energy  is  transferred  directly  from  energy-containing  to  dissipation  range. 
If  the  nonlinear  interaction  does  produce  such  a  cascade  process,  it  is 

clear  that  the  characteristic  time  associated  with  a  mode  in  the  inertial  range 

-2/3 
will  have  a  value  proportional  to  k  '  "^ ,     This  result  may  be  obtained  by 

forming  a  time  froni  the  length  k"  and  the  rms  velocity  corresponding  to  the 

portion  of  the  total  energy  which  is  contained  in  a  band  of  width  k  centered  on  k, 

According  to  the  Kolmogorov  dimensional  considerations,  it  follows  that  in  the 

inertial  range  U(k,-c)  must  be  of  the  fonn 

(9.12)  U(k,x)  =  cc^/^  k-^/3  r(e^/V/^x), 

where  c  is  a  numerical  constant,  e  is  the  total  rate  of  energy  transfer/unit  mass 
by  the  assumed  cascade  process,  and  r(  )  is  a  vmiversal  function  normalized  by 
r(0)  =  1. 

We  may  now  verify  from  (9.8)  that  (9.12)  is  consistent  with  the  existence 
of  an  energy  cascade  process.  We  consider  the  contirLbution  to  S(k^,0)  due  to 
interactions  with  all  the  modes  of  wave  number  k  >  k^  »  k, ,  assuming  both  k, 
and  kp  are  well  within  the  inertial  range.  This  is  essentially  equivalent  to 
considering  the  contribution  to  S(k^ ,0)  due  to  integration  over  the  region  k'  >  kj 
(for  which  the  k"  integration  is  confined  to  k"~k').  Noting  from  (9.12)  that 
U(k«,0)  «  U(k,,0)  if  k'  >  kg,  then  using  (9.10)  to  estimate  the  relative  magnitudes 
of  the  two  terms  of  (9.8),  and  neglecting  the  exponential  factors  in  accordance 
with  the  assiuned  smallness  of  the  dissipation,  we  may  verify  that  this  contribu- 
tion ia 

(9.13)  -  [s(k^,0)]    ^|Ec2ek^-3    j  (r(x))2dxj  I  j^  I 


•^ 


o 
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Integrating  (9.13)  over  all  modes  from  the  bottom  of  the  inertial  range  (wave 

©  ha'v 
-U/3 


number  k  ,  ),  which  we  assume  is  well  below  k- ,  we  have 


<-  1  /k  \"^'-' 

(9.II4)     -  Un  J    [S(k,0)]^   k^dk  -  cc  (^ 

min 
where  c'  is  a  new  numerical  constant. 

This  result  indicates  that  if  k^A-i  is  sufficiently  large  there  is 
negligible  direct  transfer  of  energy  from  modes  of  wavenumber  <  k,  to  modes  of 
Wavenuniber  >  k_.  Thus  the  asymptotic  validity  of  the  cascade  concept  is 
corroborated. 


Further  analysis  shows  that  c'  in  (9.1U)  is  of  order  unity.  This  indicates 

that  the  number  of  distinct  cascade  steps  in  the  entire  inertial  range  is  of  the 

order  of  magnitude  of  logT„(k _A  ,•«)*  where  k „  defines  the  top  of  the  range, 

iu    max    mxn  '  max 

The  similarity  theory  for  high  Reynolds  niimber  indicates  the  rough  relation 
ln(k:       A  •    )  '^  r  In(Re),  where  Re  is  the  fteynolds  number  based  on  macroscale 
and  rotal  rms  velocity.     We  may  estimate  crudely  that  for  Re   ~  ID     there  would 
be  only  about  four  distinct  steps. 

The  rather  ill-defined  nature  of  the  cascade  process  suggests  that  even 
in  quite  low  Reyiiolds  number  turbulence  the  modes  in  the  energy-containing  range 
interact  with  many  other  modes.     This  adds  support  to  the  validity  of  the  weak 
dependence  hypothesis  in  this  case. 

9.2     Non-zero  magnetic  field.     Extending  our  results  to  isotropic  hydroma.^netic 
turbulence,  we  find 

(9.15)  W(k,x)  +  vkVk,x)     -     Q(k,T;), 

(9.16)  U(k,^)  +  vk^U(k,t)     -     S(k,t)  +  Q(k,t)  +  G(k,t), 
where 

(9.17)  <  Wj^(k,t  +  Ti)w*(k,t)  >     -    5  PijCk)  W(k,t), 

00         00  k+k» 

(9.18)  Q(k,t)  =  nk  f    ds  )     k'dk'    r       k"dk"  r^^{jj{^'^^"-^  exp(-vk^s)W(k,s) 

XW(k;s+x)TJ(k",s+L)  -  ^^^(^l^o'^'^   exp(-vk«2s)w(k,s--c)W(k«,8) 
XU(k",s)  -  ^.^{^ji^^^^J^   exp(-vk'2g)w(k,s--o)W(k",s)U(k',s)l, 
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00  (X)  k+k ' 

(9.19)  Q(k,t)  «  nk  j     ds  f    k'dk'    f        k"  dk"  [  ^^^^fk^o)'^  exp(-vk^3)0(k,8)W(k',s^t) 

''o         '^o  |k-k'| 

XW(k",s+Ti)  -  ''%l^][^]^  exp(-vk«2)u(k,s-^)W(k',s)W(k",s)]. 
Tbe  new  angular  dependence  factors  are  given  by 

(9.20)  c(k,k',k'0     =     2k"-'-k«   z(l-y^) 
d(k,k',k")      -     2(1  +  xyz) 

h(k,k',k")     =     2k"\'(l  +  xy  +  y^z  -  z^) 
j(k,k',k")     =     2k"-'-k'(z  -  3y  -  x^z  -  2xyz^) 

U(k,-c)  has  its  previous  meaning,  and  both  S(k,x)  and  G(k,Ti)  are  given  by  the 
previous  expressions.  We  note  that  k  Q(k,0)  is  proportional  to  the  net  rate  of 
energy  transfer  to  all  magnetic  modes  of  wave  number  k  as  a  result  of  the  nonlinear 
interaction, 

10.  Discussion 

The  fundamental  condition  for  the  validity  of  the  theory  we  have  presented 
is  that  each  excited  mode  interact  with  sufficiently  many  others  that  the  effect 
of  each  elementary  interaction  can  be  considered  a  siriall  perturbation.  We  have 
pointed  out  that  this  condition  should  be  met  even  for  the  energy-containing  range 
of  quite  moderate  Reynolds  number  turbulence.  If  our  subsequent  perturbation 
hypothesis  is  justified,  as  indicated  by  the  evidence  presented  in  Section  5j 
then  the  principal  limitation  to  the  physical  applicability  of  the  theory  would 
appear  to  lie  in  the  artificiality  of  our  treatment  of  the  forces  which  supply 
energy  to  the  turbulence. 

It  would  be  more  realistic  to  replace  the  f  in  (3.7)  by  terms  giving  the 
interaction  of  the  fluctuating  flow  with,  say,  a  steady  shear  flow.  It  may  be 
seen  that  this  would  have  the  effect  of  introducing  a  linear  couplin,^;  among  the 
q^,   which  would  be  incompatible  with  the  homogeneity  conditions  (U.7).  The  precise 
nature  of  the  driving  forces  should  have  comparatively  little  effect  on  the  inertial 
and  dissipation  ranges  (except  for  the  residual  anisotropy  at  high  wave  numbers), 
but  it  seems  clear  that  a  proper  description  of  the  energy-containing  range  in 
shear  turbulence  by  our  methods  requires  relaxation  of  the  homogeneity  condition 

The  assinned  homogeneous  forces  may  be  visualized  as  arising  from  a  random  volume 
distribution  of  stirring  devices. 
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along  at  least  one  axis.  The  complications  associated  with  Lhis  generalization 
are  more  severe,  in  certain  respects,  than  those  associated  with  the  relaxation 
of  stationarity  required  for  the  treatment  of  the  decay  of  homogeneous  turbulence. 
Other  generalizations  of  interest  are  the  inclusion  of  external  magnetic  fields 
and  compressibility  effects. 

It  may  be  of  value  to  re-emphasize  the  difference  between  the 
hypothesis  of  weak  dependence,  introduced  in  Section  li,  and  the  hypothesis 
employed  by  Heisenberg^  -'  and  in  somewhat  different  form  by  Chandrasekhar '-  -', 
that  the  second  and  fourth  moments  of  the  (two-time)  distribution  of  the  spatial 
Fourier  components  of  the  fields  are  related  as  in  an  independent  distribution. 
For  brevity,  we  shall  confine  discussion  to  the  idealized  dissipationless 
system  discussed  in  Section  6.  Setting  v  ,  f  =0,  we  may  obtain  readily  from 
(3,7)  the  relation 

,00 


(10.1)    ijo)  -  21  /apY  V  [  T.....(3)ds. 


p,^„X  »PY  a^X  J  m^ 


M.V.ll.A   '  ■     •    > 


If,   now,  we  assert  Heisenberg's  hypothesis,  we  obtain 

(1D.2)         R^(0)  .  2     r  (A        )2    I     R^    (s)R^(s)ds, 

■vdiere  (3,11)  and  (U.7)  have  been  used.  It  follows  from  the  properties  of  correlation 
functions  that  the  integrals  on  the  right  side  of  (10,2)  are  positive  quantities, 
so  that  this  hypothesis  leads  here  to  the  impossible  result  ft__(C)  >  0. 

The  essential  error  in  (ID, 2)  is  the  neglect  of  the  skew  elements 
T  .„  (s)  (p,,X,p,Y  not  equal  in  pairs)  which  occur  in  (10. l),  Although  these 
terms  individually  are  very  small,  their  number  is  greater  than  that  of  the 
non-skew  elements  by  a  factor  the  order  of  the  total  number  of  degrees  of  freedom 
of  the  system.  The  neglect  of  their  cumulative  contribution  is  not  justified  on 
the  basis  of  our  weak  dependence  h3TX)thesis.  Quite  a  different  matter  is  the 
assertion  Tq  .  (s)  '^   Rgg(s)R  (s)  for  the  non-skew  element  appearing  in  ($,3); 
this  approximation  follows  legitimately  from  the  weak  dependence  hypothesis. 
In  distinction  from  the  incorrect  result  above,  our  (6,12)  yields  ^^^(C')  ■  0, 

The  inconsistencies  associated  with  the  assumption  of  normal  relation 
between  second  and  fourth  moments  are  discussed  more  fully  in  [l2j }   it  is  shown 
there  that  the  implied  violation  of  the  equations  of  motion  is  a  serious  one. 
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The  immediate  program  following  from  the  work  presented  here  involves  the 
determination  of  the  spectrum  in  wave  number  and  frequency  of  isotropic  hydro- 
magnetic  turbulence  fron  equations  of  the  type  presented  in  Section  9j   and 
examination  of  the  energy  transfer  among  the  various  kinetic  and  magnetic  modes. 
This  program  entails  validation  of  assumption  (7.5),   and  the  development  of 
systematic  approximation  schemes  for  the  solutions  of  the  difficult  nonlinear 
integro-differential  equations  of  motion  for  the  covariances. 

The  more  fundamental  program  suggested  by  the  work  reported  lies  in  the 
investigation  of  the  distribution  functional  defined  in  Sections  h  and  8,   further 
development  of  the  statistical  mechanical  method  there  outlined,   and  the  applica- 
tion of  this  method  to  other  problems  of  the  many-body  type.     Of  particular 
interest  is  the  applicability  of  the  present  approach  to  quantum  statistical 
mechanics.     The  time-dependent  Schroedinger  equation  of  a  quantum  mechanical 
system  constitutes  a  set  of  first-order  equations  of  motion  for  the  orthogonal 
components  of  the  state-vector  in  any  given  representation.     In  problems  where 
an  appropriate  perturbation  treatment  is  valid  one  may  obtain  equations, 
analogous  to  those  derived  in  Sections  6  and  7,  which  determine  the  correla- 
tion and  response  matrices  of  the  orthogonal  components, 
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Addendum  to  Section  9 

VJe  shall  present  here  the  full  set  of  equations  determining  the  im- 
pulse response  functions  and  covariance  scalars  when  assur^tion  (7»5)  is  not 
made.  Consider  first  non-magnetic  turbulence.   In  place  of  (9»7)  we  have,  in 
the  general  case, 

(A.l)         G(k,r)  =  I  g(k,s)  F(k,s+T)ds  , 

and  in  place  of  (9»8)  we  have 

(A,2)   S(k,T)   =  Y  (^  I  d\'[a(k,k',|k-k*|)g(k,s)U(k',s+-c)U(  |k-k'l,s+^  ) 

-b(k,k',|k-k'|)U(k,s-r)g(k',s)U(|k-k'|,&)j  , 

over 
where  we  also  have  re-expressed  the  integration/wave  vector  space.  To  obtain  a 

complete  set  of  equations  we  must  now  express  (7»3)  for  the  isotropic  case.  From 

(3 •I)  we  may  obtain  the  equation 

(A.3)        -G(k, -c)  +  vk^G(k,x)  =  H(k,  t  )  +  F(k,  t  )  , 

where 

(A.I4)        H(k,T)  =  Im  Jk^TI  <^j(lS**t)^j(lJ-is',  t)F^(k,t+T)>  I  . 

Evaluating  H(k, r )  by  our  perturbation  procedure,  and  following  the  argument  by 

which  (7.3)  was  obtained  from  (7.1),  we  find 

2     c^      f 
(A.5)   g(k,T)+  vk\(k,T)  =  -  ^  I  ds  I  d\'b(k,k',|k-k'|)  g(k,T-s)g(k',s)U(|k-kI,s), 

(  T  >  0)  , 

upon  making  the  transformation  to  the  Fourier  integral  representation.  Equations 
(9«2),  (A»5),  (A.l),  and  (A.?)  now  form  a  complete  set  which  determine  g(k, t  )  and 
U(k,  t)  under  the  boundary  conditions  g(k,+0)  =  1,  U(k,  cd)  =  0  • 

We  shall  denote  the  imptilse  response  function  for  magnetic  modes  by 
g^(k,  t).  Instead  of  (9«18),  (9<.19)  ,  we  have  in  the  general  case 


-Ui  - 


2   '^ 


(A.6)  Q(k,T)  =  ^  j  ds  I  d\'|d(k,k',|k-fe'|)gj^(k,s)W(k'ys+T)u(|k-k'l,s+T) 

-  h(k,k',|k-k'l)W(k,s-r)gj^(k',s)U(lk-k'|,s) 

-  j(k,k',|k-k'|)W(k,s-T)g(k',s)W(|k-k'|,s)]  ^ 
and 

(A.7)  Q(k,T)  =  |-  I  ds  j  dV[a(k,k',|k-k'|)g(k,s)W(k',s+-tr)W(lk-k'l,s+T) 

-  c(k,k',lk-k'l)U(k,s-T)g^(k',s)W(|k-k'l,s)j  , 
while  the  equations  for  g(k,  t  ) ,  g^(k,  t)  are 

(A.8)        g(k,T)  +  vk^g(k,T)  -  n(k,  t:)  +  in(k,-c)  ,      (  ^  >  O)  , 

(A.9)         ^(k,T)  +  ^^gj^(k,  t)  =  m(k,T)  ,  (  T  >  0)  , 

where 

(A.IC)  m(k,T)  '  -Y     I  ds  j  d\'c(k,k',Ik-k'|)g(k,T-s)gj^(k',s)W(Ik-k'|,s)  , 

(Aoll)  m(k,T)  -  -  Y  I  ^2  I  d^k'gjj^(k,  T-s)[h(k,k',|k-k'|gjjj(k',s)U(|k-k'|s) 

-  j(k,k',lk-k'l)g(k',s)W(|k-k'l,s)]  , 

and  n(k,  t)  is  the  right  side  of  (A,5)o 

The  eqiiations  presented  for  the  impulse  response  functions  have  a 
simpler  structure  than  the  eqtiations  for  the  covariance  scalars.  Th\as,  If 
U(k,  t)  were  a  known  function,  g(k,  r  )  could  be  computed  from  (A,^  by  straight- 
ofrward  step-by-step  integration  from  the  point  t  =  0, 
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